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Abstract. We have compared composition changes of NO, 

N0 2 , h 2 o 2 , o 3 , n 2 o, hno 3 , n 2 o 5 , HN0 4 , CIO, HOCl, 

and C10N0 2 as observed by the Michelson Interferometer 30 
for Passive Atmospheric Sounding (MIPAS) on Envisat in 
the aftermath of the “Halloween” solar proton event (SPE) 
in October/November 2003 at 25-0.01 hPa in the Northern 
hemisphere (40-90°N) and simulations performed by the fol- 
lowing atmospheric models: the Bremen 2D model (B2dM) 35 
and Bremen 3D Chemical Transport Model (B3dCTM), the 
Central Aerological Observatory (CAO) model, FinROSE, 
the Hamburg Model of the Neutral and Ionized Atmo- 
sphere (HAMMONIA), the Karlsruhe Simulation Model of 
the Middle Atmosphere (KASIMA), the ECHAM5/MESSy 40 
Atmospheric Chemistry (EMAC) model, the modeling tool 
for SOlar Climate Ozone Links studies (SOCOL and SO- 
COLi), and the Whole Atmosphere Community Climate 
Model (WACCM4). The large number of participating mod- 
els allowed for an evaluation of the overall ability of atmo- 45 
spheric models to reproduce observed atmospheric perturba- 
tions generated by SPEs, particularly with respect to NO y 
and ozone changes. We have further assessed the meteoro- 
logical conditions and their implications on the chemical re- 
sponse to the SPE in both the models and observations by 
comparing temperature and tracer (CTI 4 and CO) fields. 

Simulated SPE-induced ozone losses agree on average » 
within 5% with the observations. Simulated NO y enhance- 
ments around 1 hPa, however, are typically 30% higher than 
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indicated by the observations which can be partly attributed 
to an overestimation of simulated electron-induced ioniza- 
tion. The analysis of the observed and modeled NO y par- 
titioning in the aftermath of the SPE has demonstrated the 
need to implement additional ion chemistry (HN0 3 forma- 
tion via ion-ion recombination and water cluster ions) into 
the chemical schemes. An overestimation of observed H 2 0 2 
enhancements by all models hints at an underestimation of 
the OH/H0 2 ratio in the upper polar stratosphere during 
the SPE. The analysis of chlorine species perturbations has 
shown that the encountered differences between models and 
observations, particularly the underestimation of observed 
CIONO 2 enhancements, are related to a smaller availabil- 
ity of CIO in the polar night region already before the SPE. 
In general, the intercomparison has demonstrated that dif- 
ferences in the meteorology and/or initial state of the atmo- 
sphere in the simulations causes a relevant variability of the 
model results, even on a short timescale of only a few days. 


1 Introduction 

Energetic particle precipitation has important implications 
on atmospheric chemistry. In particular, protons and asso- 
ciated electrons, generated during solar eruptions, cause spo- 
radically in-situ production of NO x and HO x radicals in- 
volved in catalytic ozone destruction. These solar proton 
events (SPEs) thus represent an important Sun-Earth connec- 
tion which contributes to the natural ozone variability. The 
quasi-instantaneous increase of odd nitrogen and hydrogen 
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due to SPEs induces perturbations of the chemical composi- 
tion of the middle atmosphere on a short-time scale. In this 
sense, SPE-induced perturbations of the atmospheric com- 
position represent an ideal natural laboratory for studying 
so stratospheric and mesospheric chemistry (see also Jackman us 
andMcPeters, 1987). 

In recent years, there have been two large SPEs (Oc- 
tober/November 2003 and January 2005) (Jackman et al., 
2008) which have been intensively observed by several in- 
65 struments on different satellite platforms, including, for ex - 120 
ample, NOAA 16 SBUV/2 and HALOE data (Jackman et al., 
2005a, b; Randall et al., 2005); MIPAS, GOMOS and SC1A- 
MACHY on Envisat (Lopez-Puertas et al., 2005a, b; von 
Clarmann et al., 2005; Orsolini et al., 2005; Seppala et al., 
ro 2004; Rohen et al., 2005); and MLS on AURA (Verroneni25 
et al., 2006). In particular, during late October and early 
November 2003, three active solar regions produced solar 
flares and solar energetic particles of extremely large inten- 
sity, the fourth largest event observed in the past forty years 
75 (Jackman et al., 2005b, 2008). During and after this event, 00 
often called “Halloween” storm, the MIPAS instrument ob- 
served global changes (e.g. in both the Northern and South- 
ern polar regions, during day and nighttime) in the strato- 
spheric and lower mesospheric composition. This includes 
so enormous enhancements in NO x , e.g., in NO and NO 2 , and 135 
large depletions in O 3 (Lopez-Puertas et al., 2005a) as well 
as significant changes in other NO y species, such as HNO 3 , 
N 2 O 5 , CIONO 2 (Lopez-Puertas et al., 2005b), and N 2 0 
(Funke et al., 2008). In addition, there also have been ob- 
86 served changes in CIO and HOC1 as evidence of perturba -140 
tions by solar protons on the HO x and chlorine species abun- 
dances (von Clarmann et al., 2005). 

Several model studies, aiming at reproducing observed 
short- and medium-term composition changes after this par- 
90 ticular event (Jackman et al., 2008; Verronen et al., 2008; us 
Funke et al., 2008; Baumgaertner et al., 2010; Egorova et al., 
2010) and evaluating SPE-induced long-term effects (Jack- 
man et al., 2009) have been carried out in the past. 

The High-Energy Particle Precipitation in the Atmosphere 
95 (HEPPA) model vs. data intercomparison initiative has iso 
brought together scientists involved in atmospheric model- 
ing using state-of-the art general circulation models (GCMs) 
and chemistry-transport models (CTMs) on the one hand and 
scientists involved in the analysis and generation of observa- 
100 tional data on the other hand. The objective of this commu-155 
nity effort is (i) to assess the ability of state-of-the-art atmo- 
spheric models to reproduce composition changes induced 
by particle precipitation, (ii) to identify and - if possible - 
remedy deficiencies in chemical schemes, and (iii) to serve 
105 as a platform for discussion between modelers and data pro- mo 
ducers. This is achieved by a quantitative comparison of ob- 
served and modeled composition changes after particle pre- 
cipitation events, as well as by inter-comparing the simula- 
tions performed by the different models. 

110 In this study we report results from the intercompari- -65 


son of MIPAS/Envisat data obtained during 26 October - 
30 November 2003, before and after the Halloween SPE, 
at altitudes between 25-75 km (25—0.0 1 hPa) with simu- 
lations performed using the following GCMs and CTMs: 
the Bremen 2d Model (B2dM) (Sinnhuber et al., 2003b; 
Winkler et al., 2009), the Bremen 3d Chemical Transport 
Model (B2dM and B3dCTM) (Sinnhuber et al., 2003a), the 
Central Aerological Observatory (CAO) model (Krivolutsky 
and Vvushkova, 2002), FinROSE (Damski et al., 2007b), 
the Hamburg Model of the Neutral and Ionized Atmosphere 
(HAMMONIA) (Schmidt et al., 2006), the Karlsruhe Simu- 
lation Model of the Middle Atmosphere (KASIMA) (Kouker 
et al., 1999), the ECHAM5/MESSy Atmospheric Chem- 
istry (EMAC) model (Jockel et al., 2006), the modeling tool 
for SOlar Climate Ozone Links studies (SOCOL and SO- 
COLi) (Egorova et al., 2005; Schraner et al., 2008; Egorova 
et al., 2010), and the Whole Atmosphere Community Cli- 
mate Model (WACCM4) (Garcia et al., 2007). Among 
the species affected by SPEs we focus here on NO, NO 2 , 
H 2 0 2 , 0 3 , N 2 0, HNO 3 , N 2 O s , HNO 4 , CIO, HOC1, and 
C10N0 2 . For these species a significant perturbation well 
above the detection limit has been observed by MIPAS. We 
have further assessed the meteorological background condi- 
tions in both the models and the real atmosphere as observed 
by MIPAS by comparing temperature and tracer fields (CH 4 
and CO). Although SPE-induced composition changes dur- 
ing the Halloween event have been reported in both hemi- 
spheres, we restrict our analysis to the Northern hemisphere 
(NH) in the latitude range 40-90°N where most pronounced 
effects have been observed and composition changes can be 
well distinguished from the background variability. 

Apart from the initial particle forcing leading to atmo- 
spheric ionization, SPE-induced composition changes are 
controlled by several other factors such as the neutral and 
ion chemistry responsible for the repartitioning of primar- 
ily generated species, the background composition interfer- 
ing with the chemical repartitioning, and the meteorologi- 
cal/dynamical conditions. The large number of controlling 
factors and their interaction introduce a significant spread in 
the model results and make their analysis difficult. In order to 
reduce the model variability and to make differences between 
the simulations more traceable, we have simplified the inter- 
comparison setup such that a common particle-induced ion- 
ization source has been used in all models. These ionization 
rates, accounting for protons (154 eV-500 MeV) and elec- 
trons (154 eV-5 MeV) have been provided by the AIMOS 
model (Wissing and Kallenrode, 2009). Different model re- 
sponses to the particle forcing are hence reduced to differ- 
ences of the intrinsic model properties, e.g. chemical and 
dynamical schemes. A major aim of this paper is the assess- 
ment of these differences and their implications on the mod- 
els’ ability to correctly describe particle precipitation effects 
which represent an important source of natural, solar-induced 
climate variability on short and mid-term scales. Addition- 
ally, conclusions on the quality of the description of the ex- 
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Table 1. Used MIPAS data versions (indicated by the last digits of the retrieval version) for all species on a daily basis within the period 
26 October-30 November 2003. 
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ternal forcing provided by the ionization model can be drawn 
from the overall agreement of the short-time response of pri-205 
marily generated constituents (i.e. NO x ). 

The paper is organized as follows: In Section 2 we give 
an overview on MIPAS observations and data products used 
in this study, followed by Section 3 describing the ioniza- 
tion model AIMOS and Section 4 describing the participat- 210 
ing global circulation and chemistry transport models. The 
intercomparison method is described in Section 5, followed 
by the discussion of the results (Section 6). 

2 MIPAS observations 

215 

The Michelson Interferometer for Passive Atmospheric 
Sounding (MIPAS) is a mid-infrared Fourier transform limb 
emission spectrometer designed and operated for measure- 
ment of atmospheric trace species from space (Fischer et al., 
2008). It is part of the instrumentation of the European Envi -220 
ronmental Satellite (ENVISAT) which was launched into its 
sun-synchronous polar orbit of 98.55°N inclination at about 
800 km altitude on 1 March 2002. MIPAS passes the equator 
in a southerly direction at 10.00 a.m. local time 14.3 times 
a day, observing the atmosphere during day and night with 225 
global coverage from pole to pole. The instrument’s field of 
view is 30 km in horizontal and approximately 3 km in ver- 
tical direction. MIPAS operated during October/November 
2003 at full spectral resolution of 0.035 cm' 1 (unapodized) 
in terms of full width at half maximum. During this pe -230 
riod, MIPAS recorded a rear-viewing limb sequence of 1 7 
spectra each 90 s, corresponding to an along track sampling 
of approximately 500 km and providing about 1000 vertical 
profiles per day in its standard observation mode. Tangent 
heights covered the altitude range from 68 down to 6 km 235 
with tangent altitudes at 68, 60, 52, 47, and then at 3 km 
steps from 42 to 6 km. 

Trace gas profiles have been retrieved from calibrated ge- 
olocated limb emission spectra with the scientific MIPAS 
level 2 processor developed and operated by the Institute of 2 ® 
Meteorology and Climate Research (IMK) in Karlsruhe to- 
gether with the Institute de Astrofisica de Andalucia (IAA) 


in Granada. The general retrieval strategy, which is a con- 
strained multi-parameter non-linear least squares fitting of 
measured and modeled spectra, is described in detail in von 
Clarmann et al. (2003). Its extension to retrievals under con- 
sideration of non-LTE (i.e., CO, NO, and NO 2 ) is described 
in Funke et al. (2001). Non-LTE vibrational populations 
of these species are modeled with the Generic RAdiative 
traNsfer AnD non-LTE population Algorithm (GRANADA) 
(Funke et al., 2007) within each iteration of the retrieval. 

In contrast to previous work describing MIPAS obser- 
vations of composition changes during the Halloween SPE 
(Lopez-Puertas et al., 2005a, b; von Clarmann et al., 2005), 
we base our analysis here on reprocessed IMK/IAA MIPAS 
data which have substantially improved with respect to pre- 
vious data versions. These improvements include updates 
in the LIB processing (version 4.61/62 instead of 4.59) per- 
formed by ESA as well as changes in the L2 processing 
performed at IMK/IAA. The new data set also offers full 
temporal coverage over the period of interest (26 October- 
30 November 2003). In the following, we summarize the im- 
provements of the retrieval setups for each species/parameter 
and characterize the data used in our analysis in terms of es- 
timated single measurement precision and vertical resolution 
obtained from the full width at half maximum of the rows 
of the averaging kernel (AK) matrix (Rodgers, 2000). The 
AK diagonal elements are also discussed as a measure of the 
sensitivity of the retrieval at a given profile grid point to the 
’’true” profile. Values close to zero (typically <0.03) indicate 
that there is no significant sensitivity to the retrieval parame- 
ter at the corresponding altitude and hence are excluded from 
our analysis. This value may appear unreasonably small but 
since IMK/IAA retrievals are not constrained by optimal es- 
timation (Rodgers, 2000) but by a first order smoothing con- 
straint using a Tikhonov (1963) formalism, low values do 
not hint at a large a priori content of the retrieval but only at 
extensive smearing of information over altitude. A detailed 
discussion of systematic retrieval errors can be found in pre- 
vious works describing the individual constituent retrievals 
which are referenced in the following. 
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Fig. 1. Temporal evolution of MIPAS temperature, CH-i, CO, NO, N0 2 , (from top to bottom) volume mixing ratio, single measurement 
precision, vertical resolution, and AK diagonal element profiles (left to right) during 26 October - 30 November 2003 averaged over 60- 
90°N. White dashed lines indicate 30, 40, 50, 60, and 70 km geometric altitude levels. White regions reflect meaningless data (AK diagonal 
elements smaller than 0.03). 


(2003) the joint retrieval of horizontal temperature gradients. 
Differences between both versions are of minor nature and do 
not noticeably affect the data characteristics. In the period 
of interest, observed temperatures range at 60-90°N from 
around 200 K in the lower stratosphere to around 270 K at 


2.1 Temperature 


Temperature data versions used here are V30_T-9 and 
V3O-T-10 (see Table 1), both including as an extension to 
the original retrieval setup described in von Clarmann et al. 
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the stratopause (see Figure 1). The single measurement pre - 300 
cision ranges from 0.5 K to 1 .5-2.5 K above the stratopause. 
Vertical resolution is 3-4 km below 1 hPa and 5-7 km above. 
Meaningful data are obtained in the whole vertical range of 
interest (25-70 km). 

2.2 CH 4 306 

We use version V30_CH4„12 (see Table 1) which has been 
jointly retrieved with N 2 0 (V3CLN2CL12). The retrieval 
setup is similar to that described in Glatthor et al. (2005). 
The single measurement precision ranges from 1 0-20 ppbv 310 
in the upper stratosphere to 50-70 K above and below (see 
Figure 1) Vertical resolution is 3-6 km below 0.03 hPa and 
slightly higher above. Meaningful data is obtained in the 
whole vertical range of interest (25-70 km). 

315 

2.3 CO 

CO data versions used here are V30.C0.9, V3CLCOJ0, 
described in detail in Funke et al. (2009), as well as the 320 
most recent version V3CLCO-11 (see Table 1). Improve- 
ments implemented in the latter version include an extended 
set of spectral fitting windows resulting in a better precision 
and vertical resolution in the lower and middle stratosphere. 

In the period of interest, the temporal evolution of M1PAS 
CO abundances at 60-90°N indicate polar winter descent of 325 
mesospheric air masses of about 1 0 km around 1 hPa (see 
Figure 1). The single measurement precision ranges from 
20-30% above 1 hPa to 70-80% in the lower stratosphere. 
Vertical resolution is 6-12 km below 0.1 hPa. Meaningful 
data are obtained in the whole vertical range of interest (25- 
70 km). 330 

2.4 NO 

We use version V30-N0-14 (see Table 1), available for the 
whole time period. This version has substantially improved 335 
with respect to the retrieval setup described in Funke et al. 
(2005) and the data discussed in Lopez-Puertas et al. (2005a) 
by i) the use of log(vmr) instead of vmr (volume mixing ra- 
tio) in the retrieval vector, ii) a revised correction scheme 
for line of sight variations of the NO x partitioning close to 
the terminator, and iii) joint-fitted vmr horizontal gradients 340 
at constant longitudes and latitudes. NO increases of several 
100 ppbv have been observed at 60-90°N during the intense 
proton forcing during 29 October-4 November in the upper 
stratosphere around 0.2 hPa (see Figure 1). Above, NO in- 
creases were mainly produced by polar winter descent of up- 
per atmospheric air masses, resulting in vmrs up to 1 pprnv 345 
below 70 km. The single measurement precision is of the or- 
der of 10%. Vertical resolution ranges from 4 to 8 km below 
70 km. Meaningful data are obtained in the whole vertical 
range of interest (25-70 km). 


2.5 N0 2 

NO 2 data versions used here are V30_N02_11, 
V3CLN02J3, and V30.N02 I4 (see Table 1). In- 
cluding the same modifications as described above for NO, 
these versions have substantially improved with respect 
to the retrieval setup described in Funke et al. (2005) and 
the data discussed in Lopez-Puertas et al. (2005a). While 
differences between the latter two versions do not affect 
noticeably the data characteristics, a modified regularization 
scheme and terminator treatment implemented after version 
V30_N02_1 1 gave rise to non-negligible differences in 
the newer versions with respect to the previous setup. 
These differences are visible in the vertical resolution in 
the mesosphere and middle stratosphere (see Figure 1, 
third column) and go along with generally smaller vmrs 
around the terminator at 70°N around 0.1 hPa. Similar 
to NO, increases of 50 — 80 ppbv were observed during 
the proton forcing in the upper stratosphere, descending 
by approximately 10 km until the end of November. Polar 
winter descent of NO x led to mesospheric NO 2 increases 
of more then 100 ppbv, particularly in the second half of 
November. The single measurement precision is of the order 
of 5-10%. Vertical resolution ranges from 4 to 8 km below 
70km. Meaningful data are obtained in the whole vertical 
range of interest (25-70 km). 

2.6 N 2 0 

We use version V30_N20_12 (see Table 1), available for the 
whole time period. This version, which has already been 
used for the previous analysis of N 2 0 abundance changes 
during the Halloween SPE (Funke et al., 2008) and dif- 
fers from other versions by a relaxed regularization above 
approximately 40 km which allows for vertically resolving 
the upper stratospheric and mesospheric enhancements. At 
60-90°N, these enhancements of around 5-7 ppbv appeared 
around 30 October and descended during November to the 
middle stratosphere (see Figure 2). The single measurement 
precision ranges from 0.5 ppbv in the upper stratosphere to 

2 ppbv in the mesosphere. Vertical resolution is 4-6 km. 
Meaningful data are obtained in the whole vertical range of 
interest (25-70 km). 

2.7 HNO 3 

We use version V30_HN03_9 (see Table 1 ), available for the 
whole time period, and which is based on the retrieval setup 
described in Wang et al. (2007). HNO 3 increases of around 

3 ppbv up to altitudes of 0. 1 hPa during the proton forcing 
and a further buildup at slightly lower altitudes at the end 
of November are visible in Figure 2, consistent with previ- 
ous findings (Lopez-Puertas et al., 2005b). The single mea- 
surement precision ranges from 0. 1 ppbv in the middle strato- 
sphere to 0.35 ppbv around the stratopause. Vertical resolu- 
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Fig. 2, Same as Figure 1, but for NoO, HNOa, N0O5, and HNO4. 

3 so tionis 3-4 km below 12 hPa and 7-1 Okm above. Meaningful below approximately 0.3 hPa(52km). 
data are obtained below 0.1 hPa (60 km). 

365 2.9 HNO 4 

2.8 N 2 0 5 

We use version V3CLHN04J2 (see Table 1) which differs 
N 0 O 5 data versions used here are V3CLN205 9 and from the original retrieval setup described in Stiller et al. 

V3O.N2O5.10 (see Table 1), all based on the retrieval setup (2007) by the application of a weaker regularization in the 

35 s described in Mengistu Tsidu et al. (2004). Differences be- middle stratosphere, where most pronounced SPE effects are 

tween both versions are of minor nature and do not affect arc expected. Unfortunately, this version is sensitive to system- 

noticeably the data characteristics. N 2 0 5 increases related atic oscillations in the radiance baseline related to an im- 

to the proton event are mainly visible in Figure 2 in the perfect gain calibration of the instrument (see also Stiller 

second half of November around 2-0.5 hPa, consistent with et al., 2008). In consequence, retrieved HNO4 profiles are 

sec previous findings (Lopez-Puertas et al., 2005b). The single systematically biased during each gain calibration period 

measurement precision ranges from 0.05 ppbv to 0.15 ppbv 375 (typically a few days) with a randomly changing magni- 

in the middle stratosphere. Vertical resolution is 5-7 km be- tude from one calibration period to another. The variable 

low 2 hPa and 5-7 km above. Meaningful data are obtained bias is noticeable in the temporal evolution of the observed 
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Fig. 3. Same as Figure 2, but for O 3 , H 2 O 2 , CIO, HOC1, and CIONO 2 . 

HNO4 distributions at 60-90°N (see Figure 2) as sharp in- 0.1 5 ppbv (hardly visible in Figure 2). The single measure- 

creases/decreases in the upper stratosphere, coincident with ment precision ranges from 0.03 ppbv to 0.1 5 ppbv around 

380 the onsets of new gain calibration periods (i.e., 28 October, the stratopause. Vertical resolution is 6-1 Okm below 5 hPa 

10 November, and 24 November). Therefore, we restrict our and around 10 km above. Meaningful data are obtained be- 

analysis of SPE-related HNO4 increases in Section 6 to data asm low approximately 0.2 hPa (55 km), 

observed during one particular gain calibration period, 28 
October - 5 November, covering the onset of the proton forc- 
ing which led to short-term HNO4 increases of the order of 


385 
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2.10 0 3 

We use version V30-03-9 (see Table 1), available for the 
whole time period. Retrieval setup and characteristics are 
similar to those described in Steck et al. (2007), except for 
395 a slightly different selection of spectral intervals (micro- 
windows) and the inclusion of pre-fitted horizontal temper- 
ature gradients. A pronounced O3 depletion during the in- 
tense proton forcing as already reported in Lopez-Puertas 
et al. (2005a) is visible above the stratopause at 60~90°N, 

400 competing in weeks with seasonal mesospheric O3 buildup 
(see Figure 3). Also, the previously reported NO x -induced 
losses at lower altitudes are seen on a mid-term scale. The 
single measurement precision ranges from 0. 1 ppmv around 
the stratopause to 0.25 ppbv above and below. Vertical reso-440 
405 lution is 3-4 km below 1 hPa and 5-7 km above. Meaningful 
data are obtained in the whole vertical range of interest (25 — 

70 km). 

445 

2.11 H 2 0 2 

We use version V30_H202_4, available for the whole time 
410 period (see Table 1) which is based on the retrieval setup de- 
scribed in Versick (2010). H2O2 increases up to 0.1 5 ppbv 450 
have been observed at 60-90°N during the intense proton 
forcing on 29 October-^! November in the upper stratosphere 
around 0.2 hPa (see Figure 3). The single measurement pre- 
415 cision in the middle stratosphere ranges from 0. 1 to 0.2 ppbv, 
being thus of the order of the observed enhancements. In 
consequence, averaging is required for the analysis. Vertical 
resolutions larger than 10 km indicate that no relevant infor- 4 “ 
mation on the vertical distribution of the middle/upper strato- 
420 spheric enhancements can be extracted from the measure- 
ments. Meaningful data are obtained below approximately 
0.2 hPa (55 km). 

460 

2.12 CIO 

CIO data versions used here are V3O_CLO_10 and , a . 
425 V30-CLO-1 1 (see Table 1), all based on the retrieval setup 

described in Glatthor et al. (2004). Differences between both 
versions are of minor nature and do not affect noticeably the 
data characteristics. As in the case of FINO4, CIO data is af- 
fected by systematic oscillations in the radiance baseline re- 
430 lated to an imperfect gain calibration of the instrument, how- 
ever, to a lesser degree than in the case of HNO4. The sin- 47C 
gle measurement precision ranges from 0.2 ppbv in the lower 
stratosphere to 0.7 ppbv around 2 hPa, being thus higher than 
1 00% at the CIO peak height (see Figure 3). In consequence, 

435 averaging is required for its analysis. Vertical resolution is 
6-1 0 km below 2 hPa and 15-20 km above. Meaningful data 
are obtained below approximately 0.5 hPa (40 km). 


species 

energy 

pressure 

approx, altitude 



(hPa) 

(km) 

protons 

154 eV-500 MeV 

10^ 7 -100 

18-240/440 

electrons 

1 54 eV-5 MeV 

10~ 7 -1 

48-240/440 


Table 2. A1MOS particle energy ranges and the corresponding pres- 
sure and altitude levels. As the upper altitude border for protons and 
electrons lies in the thermosphere, it varies by solar activity. The 
first number indicates solar minimum conditions while the second 
number represents solar maximum. 

2.13 HOC1 

HOC1 data versions used here are V30 JIOCL.3 and 
V30_HOCL_4 (see Table 1), all based on the retrieval setup 
described in von Clarmann et al. (2006). Differences be- 
tween both versions are of minor nature and do not affect no- 
ticeably the data characteristics. flOCl increases of around 
0.3 ppbv show up in Figure 3 immediately after the main pro- 
ton forcing at the beginning of November, consistent with 
previous findings (von Clarmann et al., 2005). The single 
measurement precision ranges from 0.05 to 0. 1 ppbv around 
2hPa. Vertical resolution 8-12 km below 2hPa and coarser 
then 15 km above. Meaningful data is obtained below ap- 
proximately 0.5 hPa (40 km). 

2.14 CIONO 2 

CIONO2 data versions used here are V30-CL0N02-1 1 and 
V30-CL0N02.12 (see Table 1), all based on the retrieval 
setup described in Hopfner et al. (2007). Differences be- 
tween both versions are of minor nature and do not affect 
noticeably the data characteristics. C10N0 2 increases of 
around 0.5 ppbv are visible in Figure 3 after the main proton 
forcing above 5 hPa and last until the end of November. This 
is consistent with previous findings based on data version 
Vl_CLON02J (Lopez-Puertas et al., 2005b) in qualitative 
terms, however, the peak height of the increases is slightly 
higher (~5 km) in the newer data versions included here. 
This difference is mainly related to a change of the height- 
dependent regularization strength in order to allow for more 
sensitivity at lower and higher altitudes. The single mea- 
surement precision ranges from 0.06 to 0. 1 2 ppbv, increasing 
with altitude. Vertical resolution is 5-8 km below 2 hPa and 
12-14 km above. Meaningful data is obtained below approx- 
imately 0.5 hPa (40 km). 

3 Ionization rates 

The model intercomparison is based on ionization rates cal- 
culated with Atmospheric Ionization Module OSnabriick 
(AIMOS). The reason is to avoid different model results due 
to different ionization rates as to better understand the differ- 
ences in the dynamical and chemistry schemes of the models 
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proton and electron ionization rate at 40- 90N 



26. Oct 3. Nov 11. Nov 19. Nov 27. Nov 
2003 


Fig. 4. Temporal evolution of area- weighted averages (40-90° N) 
of AIMOS ion pair production rates for protons and electrons dur- 
ing the period of interest. Electron ionization rate below 1 hPa is 
induced by bremsstrahlung only. 


under assessment. AIMOS calculates ionization rates due 
to precipitating solar and magnetospheric particles. The al- 
titude range of calculated ionization rates is defined by the 
energy range of the particles considered, which is specific to 
480 the satellite instruments used. The data used here and their 
altitude coverage are listed in Table 2. Given by the altitude 
range of this study, the focus lies on solar particles. As parti- 
cle precipitation strongly depends on the geomagnetic field, 
the model accounts for different spatial precipitation zones. 

485 A detailed description on AIMOS can be found in Wissing 
and Kallenrode (2009). 

AIMOS is composed of two parts. One describes the spa- 
tial particle flux on top of the atmosphere while the second M5 
calculates the resulting ionization rate. Both parts will be 
490 discussed in the following. 
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Fig. 5. Altitude-latitude sections of AIMOS ion pair production 
rates for protons (top) and electrons (bottom) on 28 October 2003. 


In summary, the first part of the model describes the in- 
coming particle flux at every grid point. The spatial reso- 
lution is 96 zonal cells, divided into 48 meridional sections. 
Regions of similar particle flux are combined as, e.g., the po- 
lar cap. Given by the scaling of the mean precipitation maps, 
the temporal resolution is limited by the POES orbit and has 
been set to 2 h. 

3.2 Modeling ionization rates 


3.1 Spatial particle flux 

The particle flux on top of the atmosphere is measured by the 
TED and MEPED instruments on POES 15/16 as well as the 
SEM instrument on GOES 10. As all particle measurements 
495 are in-situ, the main challenge is to derive a global cover- 
age at any time. Inside an empirically determined polar cap 
where particle precipitation is homogeneous, the high ener- 
getic particle flux from GOES and the mean flux values from 
polar cap crossings of the POES satellites are used. Outside 
see the polar cap, particle precipitation depends on geomagnetic 
latitude, geomagnetic activity and local time. Therefore, 
mean precipitation maps for the POES TED and MEPED 
channels, based on a 4 year data set, have been produced, 
sorted by the geomagnetic K p -index and local time. These 
505 mean precipitation maps represent the spatial distribution, in- 
cluding, e.g., the movement of the auroral oval. According to 
the recent K p -level, the mean precipitation maps are selected 
and scaled to recent POES particle flux. 


The second part of AIMOS is the atmospheric particle de- 
tector model, which simulates particle interactions based on 
the GEANT4-Simulation Toolkit (Agostinelli et al., 2003). 
520 GEANT4 provides Monte-Carlo based algorithms to model 
energy deposition/ionization of protons and electrons. The 
atmospheric detector model is divided into 67 logarithmi- 
cally equidistant pressure levels, ranging from sea-level to 
1.7 x 10~" 5 Pa. Since the atmospheric parameters (density, 
525 altitude, composition and temperature) depend on latitude, 
season and solar activity, model versions for 80°N, 60°N, 
60°S and 80°S, 3 different FI 0.7 flux values and 4 differ- 
ent months are used. These parameters are adopted from the 
HAMMONIA (Schmidt et al., 2006) and MS1S (Picone et al., 
530 2002) models. The ionization rates for mono-energetic and 
isotropic particle ensembles are determined. As a final step, 
the mono-energetic ionization rates are combined with mul- 
tiple power-law fits of the particle flux at various regions. 

Figure 4 shows the temporal evolution of the resulting ion 
535 pair production rates averaged over 40-90°N during the pe- 
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riod of interest. The latitudinal distribution of the ionization 
rates of protons and electrons, respectively, is shown in Fig. 5 
for the 28 October 2003. 

The ionization rates should provide a similar forcing forsw 
all models, therefore the original data set has been adopted 
to every model grid. The data set and the adoption routine for 
a used specific grid is available at http://aimos.physik.uos.de. 

595 

4 Description of participating models 

4.1 Bremen 2d model (B2dM) 

The Bremen two-dimensional model is based on the two- 
dimensional transport, chemistry and radiation model for- 600 
merly described in Sinnhubcr et al. (2009) and Chipperfield 
and Feng (2003). It uses the dynamical core of the so-called 
two-and a half-dimensional” model THIN AIR (Kinners- 
ley, 1998), which calculates temperature, pressure, and hor- 
izontal transport on isentropic surfaces, interactive with the 
model chemistry. The model covers the altitude range fromeos 
the surface to 1 00 km in 29 isentropic surfaces, providing a 
vertical spacing of about 3.5 km. The horizontal resolution 
is about 9.5 degree. Stratospheric dynamics are forced by the 
amplitudes of waves 1 to 3 of the Montgomery potential from 
meteorological analyses with a repeating annual cycle for the 6 , 0 
period of May 1980 to April 1981. There is no quasi bien- 
nial oscillation (QBO) in the model, i.e., the modeled tropical 
stratospheric wind is always in a weak easterly state. In this 
sense, the Bremen 2d model is a two-dimensional chemistry- 
climate model which is forced to repeat a very similar see- 61S 
nario by the repeating annual cycle of the Montgomery po- 
tential. 

The chemistry is based on the SLIMCAT chemistry (Chip- 
perfield and Jones, 1999), but adapted for the use in the 
mesosphere in several ways: (1) above 50km, no family ap- 620 
proach is used; (2) H 2 O and CO 2 are treated as short-lived 
species explicitly, and H 2 is varied as well, to provide a re- 
alistic description of mesospheric HO* and CO. (3) NO x 
and HO x production by atmospheric ionization is parameter- 
ized based on Porter et al. ( 1 976) and Solomon et al. (1981), 625 
i.e., 1.25 NO x are produced, of which 45% are produced as 
N, and 55% as NO, and up to 2 HO x are produced per ion 
pair depending on pressure and ionization rate, equally dis- 
tributed to H and OII. Ionization due to Galactic Cosmic 
Rays in the stratosphere has been included based on Heaps 630 
(1978); the additional ionization due to solar and magneto- 
spheric particles is considered by introducing atmospheric 
ionization rates of protons and electrons provided by the 
AIMOS model (see section 3). 

All reaction and photolysis rates are taken from Sander 
et al. (2006). The Bremen 2d model has been used in the 
past to investigate the impact of large solar proton events as 
on the composition of the middle atmosphere (Sinnhuber 
et al, 2003b; Rohen et al., 2005; Winkler et al, 2008). For 


the HEPPA intercomparison, the two-dimensional model has 
been combined with a one-dimensional model sharing the 
same description of chemistry in the following way: 

25 model runs with the two-dimensional model are carried 
out at different longitudes, to take into account the tilt of the 
geomagnetic poles. For every MIPAS measurement used in 
the intercomparison, a one-dimensional model run is started 
initialized with output of the 2-dimensional model runs in- 
terpolated to the geo-location of the measurement, at local 
noon of the day before the measurement took place. One- 
dimensional model runs are then carried out until the time of 
the measurement, providing model output at the exact time 
and geo-location of the measurement. 

4.2 Bremen 3d Chemistry and Transport model 
(B3dCTM) 

The Bremen three-dimensional Chemistry and Transport 
Model is a combination of the Bremen transport model 
(Sinnhuber et al., 2003a) and the chemistry code of the Bre- 
men 2D model (Sinnhuber et al., 2003b; Winkler et al., 
2008), which is based on the SLIMCAT model (Chipperfield 
and Jones, 1999). 

The model has 28 isentropic levels ranging from 330 to 
3402K (approx. 10-60km) and has a horizontal resolution 
of 3.75° x 2.5°. Output is provided hourly. The vertical 
transport across the isentropes is calculated through diabatic 
heating and cooling rates. These rates are calculated using 
the radiation scheme MIDRAD (Shine, 1987). The horizon- 
tal transport is driven by external wind-fields. Advection is 
calculated by using the second order moments scheme by 
Prather ( 1 986). Meteorological data, such as horizontal wind 
speeds and temperatures, are taken from ECMWF ERA In- 
terim (Simmons et al., 2006). 

The model calculates the behavior of 58 chemical species, 
using a family approach for short-lived species (HO x , NO x , 
O x , CLO x , BrO x and CHO x ). It includes about 180 gas 
phase, photochemical, and heterogeneous reactions and uses 
the recent set of recommendations for kinetic and photo- 
chemical data established by the Jet Propulsion Laboratory 
(Sander et al., 2006). 

To account for ion chemistry reactions within the neutral 
code, the production of NO x and HO x is parameterized as 
suggested by Porter et al. (1976) and Solomon and Crutzen 
(1981). Hence 1.25 N atoms and about 2 HO x are pro- 
duced per ion pair. Atmospheric ionization due to solar and 
magnetospheric particles is considered by introducing atmo- 
spheric ionization rates of protons and electrons provided by 
the AIMOS model (see section 3). 

4.3 CAO 

The Central Aerological Observatory (CAO) model consists 
of a CTM and a 3D dynamical core with a horizontal 
resolution of 10°xl0° and vertical resolution of 2 km. The 
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Table 3. Summarized description of the models involved in this study. 


Model 

type 

vertical 

domain 

vert, range 
(km) 

horizontal 

resolution 

vert. res. 
(km) 

meteorological data 
nudging 

kinetic 

data 

B2dM 

CTM/GCM (2D) 

isentropic 

-0-100 

9.5° 

-3.5 

GCM 

Sander et al. (2006) 

B3dCTM 

CTM 

isentropic 

-10-60 

3.75° x 2.5° 

-3.5 

ECMWF ERA Interim 

Sander et al. (2006) 

CAO 

CTM/GCM 

altitude 

0-90 

10“x 10° 

2 

GCM 

Sander et al. (2003) 

FinROSE 

CTM 

hybrid 

-0-65 

10° x5° 

—2 

ECMWF ERA Interim 

Sander et al. (2006) 

HAMMONIA 

CCM 

pressure 

-0-200 

4° x4° 

-3.5 

ECMWF below 179hPa 

Sander et al. (2006) 

KASIMA 

CCM 

pressure 

-7-120 

5.6° x 5.6° 

0.75-3.8 

ECMWF below 1 hPa 

Sander et al. (2003) 

EMAC 

CCM 

hybrid 

-0-80 

2.8° x2. 8° 

-1 

ECMWF below 200 hPa 

Sander et al. (2006) 

SOCOL(i) 

CCM 

hybrid 

-0-80 

3.75°x3.75° 

-2 

- 

Sander et al. (2000) 

WACCM 

CCM 

pressure 

-0-135 

l.9°x2.5° 

-1.5 

MERRA below 50 km 

Sander et al. (2006) 


CTM calculates the concentrations of 30 minor components, 
involved in 70 chemical and 35 photochemical reactions, in 
640 the range 0-90 km. Output is provided hourly. The reaction 
rate constants, absorption cross-sections, solar radiation in- 
tensity, and quantum outputs were assigned in the tabulated 
form according to Sander et al. (2003). The annual and daily 
variations of the solar zenith angle at a given point and its 
645 dependence on the height above the Earths surface were 
taken into account. For zenith angle higher than 75°, Chap- 
mans functions have been used in accordance with Swinder 
and Gardner (1967). Photolysis rates have been recalculated 
every hour during the integration of the model. A family ap- 
650 proach (Turco and Whitten, 1974) has been used for solving 
the chemical equations, including O x (0.j+0+0( 1 D)), NO y 
(N+ NO+NO 2 +NO 3 + 2 N 2 O 5 +HNO 3 +HNO 4 +CIONO 2 ), 
Cl y (Cl+ClO+OClO+ClOO+HOCl+HCl), and HO x 
(H+OH+ HO 2 + 2 H 2 O 2 ). Other long-lived species (N 2 O, 
ess CC1 4 , CFCI 3 , CF 2 C1 2 , CH 3 CI 3 , CH 4 , H 2 0, H 2 and C0 2 ) 
were included also in simulations. The CAO model applies 
additionally to electron and proton-induced ionization also 
ionization rates caused by alpha-particles provided by the' 
AIMOS model (Wissing and Kallenrode, 2009). 

sen The vertical profiles of molecular oxygen and air den- 
sity were fixed during photochemical calculations. Hetero- 
geneous removal for H 2 0 2 , HNO 3 , HC1, and HNO 4 was- 
included in the troposphere. We assumed also fixed mix- 
ing ratios for long-lived and chemical families components 
665 at lower and upper boundaries during calculations in order to 
formulate the boundary conditions. Corresponding mixing 
ratio values were taken from Park et al. (1999). An accu-, 
rate, non-diffuse method for three-dimensional advection of 
trace species suggested by Prather (1986) was used to solve 
67o the continuity equation for each transported species (families 
and long-lived species). The chemical constituents were ini- 
tialized with profiles obtained from a one-dimensional model 
(Krivolutsky et al., 2001). Wind components used for trans- 
port by advection were obtained from the simulation with? 
67s the 3D dynamical model for each day of the year. Daily av- 
eraged global zonal, meridional, vertical wind components. 


and temperature have then be used in the CTM runs. More 
details concerning the chemical 3D model can be found in 
Krivolutsky and Vyushkova (2002). This model was used 
to study the response in composition and dynamics after the 
July 2000 SPE (Krivolutsky et al., 2006). The CAO simu- 
lation included in this study covers the period 26 October- 
4 November. 

4.4 FinROSE 

FinROSE is a global 3-D chemistry transport model (Damski 
et al., 2007a). The model dynamics are from external sources 
except the vertical wind, which is calculated inside the model 
using the continuity equation. In this study FinROSE has 35 
vertical levels (0-65 km), a horizontal resolution of 10° x5° 
and uses ECMWF Interim analyses (Simmons et al., 2006) 
for dynamics. Output is provided every 3 hours. The model 
produces distributions of 40 species and includes about 120 
homogeneous reactions and 30 photodissociation processes. 
Chemical kinetic data, reaction rate coefficients and absorp- 
tion cross-sections are taken from look-up-tables based on 
the Jet Propulsion Laboratory compilation by Sander et al. 
(2006). Photodissociation frequencies are calculated us- 
ing a radiative transfer model (Kylling et al., 1997). The 
model also includes formation and sedimentation of polar 
stratospheric clouds (PSCs) and reactions on PSCs. Tro- 
pospheric abundances are given as boundary conditions and 
long-lived trace gases are relaxed towards long time trends. 
The FinROSE model applies additionally to electron and 
proton-induced ionization also ionization rates caused by 
alpha-particles provided by the AIMOS model (Wissing and 
Kallenrode, 2009). 

4.5 HAMMONIA 

The Hamburg Model of the Neutral and Ionized Atmo- 
sphere (HAMMONIA) is a 3-dimensional GCM and chem- 
istry model covering an altitude range from the surface up 
to 1.7xlO~ 7 hPa. A detailed description of the model is 
given by Schmidt et al. (2006). Simulations of particle pre- 
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cipitation effects use a modified version of HAMMONIA. 

It treats 54 photochemical, 139 bi- and termolecular, 5 ion- 
715 electron recombination, and 12 ion-neutral reactions involv- 
ing 50 neutral and 6 charged (0 + , Oj , N + , Nj, NO + , e ) 
components. Neutral bi- and termolecular reactions and the 
corresponding rate coefficients are taken from Sander et al. 
(2006). Photochemistry involves 7 ionizing and dissociat- 770 
720 ing reactions through solar irradiance of wavelengths shorter 
than Lyman-alpha using a parametrization of Solomon and 
Qian (2005) and observed solar spectral irradiance. Addi- 
tionally, 6 ionizing, dissociating, and exciting reactions rep- 
resent the direct influence of precipitating primary and sec- 775 
725 ondary particles on thermospheric chemistry. Corresponding 
reaction rates are calculated using the particle induced ion- 
ization rates and branching ratios given by Roble and Ridley 
(1987) and Rusch etal. (1981). Below 10 _i hPa, particle im- 
pact on chemistry is represented by the production of N( 2 D), 7 m 
730 N( 4 S) and FIO x . Here, HAMMONIA uses parametrizations 
of Jackman et al. (2005a) based on formulations of Porter 
et al. (1976) and Solomon et al. (1981). The simulations use 
67 pressure levels and a horizontal resolution of T31 (96 lon- 
gitudes x 48 latitudes). Output is provided every 2 hours. 7 «s 
735 Up to 179hPa, the model is relaxed to ECMWF analyzed 
temperature, divergence, vorticity, and surface pressure. 

4.6 KASIMA 


proton hydrates (de Zafra and Smyshlyaev, 2001) has been 
modied to be dependent on actual ionization rates. 

4.7 EMAC 

The ECHAM5/MESSy Atmospheric Chemistry (EMAC) 
model is a numerical chemistry and climate simulation sys- 
tem that includes sub-models describing tropospheric and 
middle atmosphere processes and their interaction with 
oceans, land and human influences (Jockel et al., 2006). 
It uses the Modular Earth Submodel System (MESSy, see 
Jockel et al., 2005) to link multi-institutional computer 
codes. The core atmospheric model is the 5th genera- 
tion European Centre Hamburg general circulation model 
(ECHAM5, Roeckner et al., 2006). Here, EMAC (ECHAM5 
version 5.3.02, MESSy version 1.8 + ) was applied in the 
T42L90MA-resolution, i.e. with a spherical truncation of 
T42 (corresponding to a quadratic Gaussian grid of approxi- 
mately 2.8 by 2.8° in latitude and longitude) with 90 vertical 
hybrid pressure levels up to 0.01 hPa. The frequency of out- 
put is every 2 hours. The model is weakly nudged at 200 
- 700 hPa to ECMWF reanalysis data. The chemistry sub- 
model includes 104 chemical species and 250 homogeneous 
and heterogeneous reactions based on Sander et al. (2006). 
For more details on the setup used here refer to Baumgaert- 
ner et al. ( 2010 ). 


The KASIMA model is a 3D mechanistic model of the mid-;® 
die atmosphere including full middle atmosphere chemistry 
740 (Kouker et al., 1999). The model can be coupled to specific 
meteorological situations by using analyzed lower boundary 
conditions and nudging terms for vorticity, divergence and 
temperature. Flere we use the version as described by Red- 
dmann et al. (2010). It has a horizontal resolution of about 795 
745 5.6° x 5.6° with 63 pressure levels between 7 and 1 20 km 
and a vertical resolution in the lower stratosphere of 750 m, 
gradually increasing to 3.8 km at the upper boundary. The 
frequency of output is every 6 hours. The model is nudged 
to ECMWF analyses below 1 hPa. A numerical time step of boo 
750 1 2 min was used in the experiments. The chemistry is calcu- 

lated up to 90 km, above which only transport is applied. The 
chemistry uses JPL 2002 data (Sander etal., 2003) and is cal- 
culated up to 90 km, above which only transport is applied. 
The chemical fields of long-lived tracers have been initial- aos 
755 ized from a multi-annual run starting in the year 1960. For 
the HEPPA experiments, the transport scheme has been re- 
vised to allow transport of the members of chemical families 
NO x and HO x individually in the mesosphere. In addition, 
the ozone heating rate is calculated interactively. The rateaio 
780 constants of the gas phase and heterogeneous reactions are 
taken from Sander et al. (2003). For the production of HO x 
the parameterization of Solomon et al. (1981) is used, for the 
production of NO x , 0.7 NO molecules are produced per ion 
pair and 0.55 N atoms in ground state, including reactionsais 
765 of N + O2, N + NO, N-F N0 2 . The HNO 3 production from 


4.8 SOCOL and SOCOLi 

SOCOL (modeling tool for SOlar Climate Ozone Links stud- 
ies) is a combination of the GCM MA-ECHAM4 and the 
chemistry-transport model (CTM) MEZON. It is a spectral 
model with T30 horizontal truncation resulting in a grid spac- 
ing of about 3.75; in the vertical direction the model has 39 
levels in a hybrid sigma-pressure coordinate system spanning 
the model atmosphere from the surface to 0.01 hPa. Time 
step for dynamical and physical processes is 15 min and 2 
hours for radiative transfer calculations and chemical reac- 
tions. The original version of the chemistry-climate model 
SOCOL was described by Egorova et al. (2005), and up- 
dated version in Schraner et al. (2008). For the HEPPA com- 
parison, two versions of SOCOL have been used. One is 
with parameterized production of odd nitrogen, where for the 
NO x sources, the fact that 1 .25 NO molecules were produced 
were taken into account (Porter et al., 1976), for the HO x 
sources, the table given by Solomon et al. (1981) has been 
used. The other version (SOCOLi ) includes the chemistry 
of ionized species. SOOCLi is described in Egorova et al. 
(2010). As sources for ionization the model uses galactic 
cosmic rays (Heaps, 1978), energetic electron precipitation, 
solar proton events and observed solar irradiance. SOCOLi 
takes into account 580 reactions involving 43 neutral of the 
oxygen, hydrogen, nitrogen, carbon, chlorine and bromine 
groups, electrons, 3 1 positive and 1 7 negative charge species 
including clusters of 0‘2 + , H ‘ and NQA The rate constants 
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Fig. 6. Effect of application of averaging kernels (Alts) to the model data on the example of MIPAS and WACCM4 HOC1 zonal mean 
distribution (40-90°N) averaged over the period 29 October to 4 November 2003. 


of the gas phase and heterogeneous reactions are taken from sso 
Sander et al. (2000). SOCOL and SOCOLi models apply ad- 
ditionally to electron and proton-induced ionization also ion- 
ization rates caused by alpha-particles provided by AIMOS 
model (Wissing and Kallenrode, 2009). This choice is based 
on the assumption that AIMOS model describe all physical ess 
processes relevant to particle precipitation during the event. 
Output is provided at the local time and location of the MI- 
PAS overpass. 

4.9 WACOM 860 

The fourth version of the Whole Atmosphere Commu- 
nity Climate Model (WACCM4) is part of the Com- 
munity Earth System Model (http://www.cesm.ucar.edu/). 

It is a coupled chemistry climate model with hori- 
zontal resolution of 1.9° latitude by 2.5° longitude. 
For this study WACCM4 has 88 vertical levels and 
is forced with meteorological fields from the Mod- 
em Era Retrospective-analysis for Research and Applica-^ 
tions (MERRA, http://gmao.gsfc.nasa.gov/research/merra/). 
MERRA is a NASA reanalysis for the satellite era using the 
Goddard Earth Observing System Data Assimilation System 
Version 5 (Rienecker et al., 2008). The forcing is achieved by 
relaxing horizontal winds and temperatures with a time con- 6 __ 
stant of approximately 50 hours from the surface to 40 km. 
Above that level the forcing is reduced linearly, so that the 
model is free-running between 50 km and the model top at 
approximately 135 km (4.5x 10“ 6 hPa). Heating rates and 
photolysis are calculated using observed daily solar spec- 
tral irradiance and geomagnetic activity effects in the auroral 
region are parameterized in terms of the Kp index (Marsh 
et al., 2007). A description of simulations of the effects em 
of solar proton events using an earlier free-running version 
of WACCM and comparison with measurements is given in 


Jackman et al. (2008, 2009). The standard WACCM chem- 
istry is described and evaluated extensively in WMO (2010). 
Reaction rates are from Sander et al. (2006). For these sim- 
ulations we have modified the N+N02 reaction to include 
two additional pathways as described in Funke et al. (2008). 
It should be noted that both WACCM and HAMMONIA use 
the same chemical solver based on the MOZART3 chem- 
istry Kinnison et al. (2007), include the same set of ion- 
ized species, and use the parameterized EUV ionization rates 
from Solomon and Qian (2005). For these simulations the 
latter parameterization has been extend to include the pho- 
toionization of COa in the EUV. Proton and electron ion- 
ization rates, used in the nominal simulation, are taken from 
AIMOS, however above 5xl0- 4 hPa (MOO km) ionization 
from electrons is instead calculated by the WACCM param- 
eterized aurora. An additional simulation using proton ion- 
ization, only, has also been performed (in the following de- 
noted as WACCMp). The HO x production per ion pair is in- 
cluded in WACCM using a lookup table from Jackman et al. 
(2005b, Table 1), which is based on the work of Solomon et 
al. (1981). It is assumed that 1.25 N atoms are produced per 
ion pair and divide the N atom production between ground 
state, N( 4 S), at 0.55 per ion pair and excited state, N( 2 D), at 
0.7 per ion pair (Jackman et al., 2005b; Porter et al., 1976). 
In this study, WACCM constituent and temperature profiles 
were saved at the model grid point and time-step (model 
time-step is 30 minutes) closest to each of the MIPAS ob- 
servation locations. 

5 Intercomparison method 

In order to reduce errors related to the different sampling 
of the MIPAS observations and gridded model data (i.e., 
B3dCTM, CAO, FinROSE, HAMMONIA, KASIMA, and 
EMAC) , we have linearly interpolated the model results to 
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Fig. 7. Observed and modeled temperature zonal mean distributions at 40-90°N averaged over the period 29 October to 4 November 2003. 


the MIPAS measurement locations and times, as well as to 
the corresponding pressure levels of the vertical retrieval gridsis 
ass of the species under consideration. This approach has the fur- 
ther advantage that diurnal variations of particular species are 
implicitly taken into account. Comparison of MIPAS mea- 
surements and model results requires the transformation of 
modeled profiles to MIPAS altitude resolution. Based on the 
890 formalism by Rodgers (2000), we calculate the model pro- 
files adjusted to MIPAS resolution x a dj as 

920 

Xailj — Ax mod + (J A)X a , ( 1 ) 

where A is the MIPAS averaging kernel matrix, x mo d is the 
original model profile, I is unity, and x n is the a priori in- 
ass formation used in the MIPAS retrievals. Assuming that the 925 
altitude resolution of the models is much finer than that of 
the MIPAS retrievals, the comparison of x a dj and MIPAS 
measurements is not affected by any smoothing error. 

This procedure has been applied to each model result re- 
900 sampled at the corresponding measurement location. As 930 
an example. Figure 6 compares HOC1 zonal mean distribu- 
tions at 40-~90°N, averaged over the period 29 October to 
4 November 2003, as observed by MIPAS and as modeled 
by WACCM with and without application of averaging ker- 
305 nels. In the latter case, the vertical distribution is broader 935 
and slightly shifted towards lower altitudes, similar to the re- 
trieved MIPAS profiles. Also, the absolute vmr peak values 
are smaller than without application of the averaging kernel. 

It should be noted that the apparent better agreement of the 
910 maximum vmr values between MIPAS and the unconvolved <mo 
WACCM simulations are related to the fact that background 
HOC1 vmrs are underestimated in the model. The relative 
vmr increase related to the SPE is in better agreement when 


comparing observed and convolved model data (see also Sec- 
tion 6 ). 


6 Results and discussion 

6.1 Meteorological background conditions 

Meteorological background conditions, particularly the ther- 
mal structure and the prevailing dynamics, can have an im- 
portant impact on the magnitude and spatial distribution of 
SPE-induced composition changes. Temperature differences 
between models and observed data have, on the one hand, 
a significant impact on SPE-related chemistry due to in- 
volvement of highly temperature-dependent reactions (i.e., 
N+0 2 or NO 2 +O 3 ). On the other hand, meridional trans- 
port and mixing, depending largely on the development stage 
of the early winter polar vortex, control the redistribution 
of air masses between polar night and illuminated regions 
and hence, the efficiency of photochemical losses. Also, the 
strength of polar winter descent plays an important role in 
the vertical redistribution of some species on the lime scale 
of the intercomparison period. 

Figure 7 shows the MIPAS temperature zonal mean dis- 
tribution at 40-90°N averaged over the period of the main 
proton forcing, 29 October to 4 November 2003, and the 
corresponding differences between the models and the ob- 
servations. It is evident that models which are driven or 
strongly forced by assimilated meteorological data up to 
the upper stratosphere (i.e., B3dCTM, FinROSE, KASIMA, 
and WACCM) reproduce reasonably well the observed tem- 
peratures below approximately 1 hPa. On the other hand, 
free-running models (B2dM, CAO and SOCOLi) and those 
which are nudged to meteorological in the troposphere. 
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Fig. 8. Temporal evolution of MIPAS temperature and model - MIPAS differences averaged over 70-90°N. 
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Fig. 9. Relative meridional CH 4 anomaly in MIPAS observations and model simulations at 40-90°N averaged over the whole time period. 


only (EMAC and HAMMONIA) tend to overestimate the 
observations inside the polar vortex by more than 15 K. 
around approximately 1 hPa or slightly below. Slightly too 
high stratopause temperatures are found in EMAC, CAO, 
B3dCTM, and FinROSE simulations. In the polar meso-^ 
sphere, temperatures are underestimated by HAMMONIA 
(up to 15 K) and SOCOLi (more than 25K). The temporal 
evolution of observed polar temperatures (70-90°N) and the 
corresponding differences between model and observations 
are shown in Figure 8. No significant trend in both, observa- gt6 
tions and model data, can be observed during the period of 
interested, while short-term temperature fluctuations of more 


than 10 K compared to the observations, most likely related 
to differences in the planetary wave activity, show up particu- 
larly in the case of the free-running or weakly nudged models 
(B2dM, CAO, HAMMONIA, and SOCOLi). 

Differences in the magnitude of meridional redistribu- 
tion between models and observations have been assessed 
by comparing CH4 zonal mean distributions provided by 
all models except CAO. Since the global stratospheric CH4 
abundances differ noticeably among the models, we used the 
relative meridional CH,| anomaly as indicator for meridional 
redistribution rather than absolute vmrs. Figure 9 shows the 
observed and modeled meridional CH4 anomalies at 40- 
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Fig. 1 1. Temporal evolution of CO abundances in M1PAS observations and model simulations averaged over 70-90°N. The 0.1, 0.2, 0.5, 1. 
2, 5, 10, and 20 ppmv contours are shown by solid lines. 


90° N averaged over the whole period. A pronounced gra- 
dient in the observed anomalies around 60°N indicates the 
870 early winter vortex boundary. In general, the vortex bound - bm 
ary position is well reproduced by all models, although there 
are significant differences in the overall CH.j change from 
mid-latitudes to the pole between the models. Strongest lat- 
itudinal gradients (i.e., weakest redistribution) were found 
975 in the KASIMA simulations, while smallest gradients (i.e., ass 
strongest redistribution) are visible in HAMMONIA. The 
reason for the underestimation of meridional redistribution 


in KASIMA is not fully understood, particularly because 
other ECMWF-driven models have simulated considerably 
stronger mixing. The too strong mixing in HAMMONIA 
is most probably related to wave-1 activity, being present in 
the whole time period. The vertical distribution of the ob- 
served CH 4 meridional anomaly shows a broadening in the 
stratopause region (1-0. 1 hPa), indicating a weakened trans- 
port barrier at the vortex top. This behavior is reproduced by 
the models in general although there are differences with re- 
spect to the altitude and magnitude of the broadening region. 
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In B2dM, it is shifted slightly upwards while the opposite is 
observed in EMAC, SOCOL, and SOCOLi simulations. In 
990 these latter models, meridional redistribution seems also to 
be slightly overestimated around the stratopause. It should 
be noted that our analysis of CH 4 meridional anomalies does 
not allow to distinguish between meridional redistribution by 
eddy diffusion and large-scale transport by planetary waves, 

995 the latter being of higher importance for the redistribution of 
air masses between polar night and illuminated regions. 

The variability of the polar vortex strength has been as- 
sessed by comparing the temporal evolution of the relative 
change of CH 4 abundances with respect to 26 October aver- 
1000 aged over 70-90°N (see Fig. 10). The observed evolution 
indicates a vortex intensification and descent in the lower 
and middle stratosphere while a CH 4 increase above 0.3 hPa, 
particularly during the proton forcing at the beginning of 
November, hints at an increase of meridional mixing in the 
1005 mesosphere. This general behavior is qualitatively repro- 
duced by the models, although important differences with 
respect to the vertical structure and magnitude exist. These 
differences have to be taken into account when analyzing 
the temporal evolution of SPE-induced composition changes 
1010 (see next sections). 

Carbon monoxide is an ideal tracer for upper stratospheric 
and mesospheric dynamics. Particularly, it allows to identify 
air masses which have descended from the upper mesosphere 
and contain enhanced NO x related to energetic electron pro- - 
1015 cipitation (EEP). Since polar winter descent of NO x gener- 
ated by EEP prior to the SPE event, is not resolved by all 
models and since we focus here on SPE-related effects, ob- 
served NO x enhancements due to descending upper meso- 
spheric air masses perturb our analysis and should hence baoso 
1020 excluded. MIPAS CO observations provide an excellent cri- 
terion for identification of EEP-related enhancements (see 
Section 6.2). 

CO distributions also allows for the characterization of de- 
scent and vortex perturbations by large-scale wave activityioss 
1025 and isentropic mixing across the vortex boundary in the up- 
per stratosphere and mesosphere. Figure 1 1 shows the ob- 
served and modeled temporal evolution of CO abundances 
average over 70-90°N. In general, the continuous decrease 
in altitude of CO vmr isolines in the upper stratosphere,ioeo 
1030 related to polar winter descent, agrees well in models and 
data. Around 1 hPa, polar air masses descended approxi- 
mately 5 km in both models and observations during the time 
period under investigation. 

A higher variability is found in the mesosphere. Ob- 
1035 served CO abundances decreases around the beginning ofc® 
November, at the same time when CH 4 increased signifi- 
cantly (see Fig. 10). A pronounced CIO increase occurred 
around 20 November, hinting at enhanced descent and vortex 
intensification. Modeled CO distributions show a different 
too temporal evolution in the mesosphere, although some simi- 
larities can be found. For instance, EMAC, KASIMA andnm 
WACCM reproduce the CO increase in late November, how- 



Fig. 12. Area- weighted averages (40-90° N) of observed and mod- 
eled NO y enhancements during 30 October-1 November with re- 
spect to 26 October (left) and relative deviations of modeled aver- 
ages from the MIPAS observations (right). Thick solid and dashed 
lines represent model multi-model mean average and MIPAS obser- 
vations, respectively. WACCMp denotes the WACCM simulation 
including proton ionization, only (excluded from the multi-model 
mean). 


ever, with a smaller magnitude and slightly shifted in time. A 
CO decrease around 0. 1 hPa at the beginning of November, 
as observed by MIPAS, is also visible in SOCOL, SOCOLi, 
and - to a lesser extent - in WACCM simulations. 

In contrast to the observations, these modeled decreases 
occur nearly instantaneously on 3 1 October, suggesting that 
the simulated CO changes might be related to the proton 
event rather than dynamical modulations. Indeed, CO is re- 
moved by the reaction with OH, which is strongly enhanced 
during the SPE at nighttime. The isolation of a possible SPE- 
induced chemical CO loss from dynamical effects is diffi- 
cult in both, observations and simulations. Nevertheless, we 
have analyzed the observed CO abundances at fixed CII 4 
levels in the vertical range of 0.2-0.05 hPa in order to exclude 
CO variations related to isentropic mixing or meridional re- 
distribution. CO abundances observed simultaneously with 
CH 4 vmrs of less than 40 ppbv decreased by approximately 
1 ppmv from 29 October to 1 November, thus suggesting a 
chemical removal of the order of 10 % which could be related 
to enhanced OH. The CO decreases found in the WACCM 
simulations have a similar magnitude, while SOCOL and 
SOCOLi simulations show a CO decrease around 30%. 

6.2 Enhancements of NO y and N 2 0 

The most important impact of proton precipitation on the 
middle atmosphere is the immediate formation of NO x (= 
NO + NO 2 ) via dissociation of molecular nitrogen by ion- 
ization and subsequent recombination with oxygen. Due to 
its relatively long chemical lifetime in the stratosphere, SPE- 
induced NO x enhancements have a strong potential to de- 
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plete ozone on a tnid- to long-term scale via catalytic cycles. 

A fraction of excess NO x produced by proton forcing is sub- 
sequently buffered into NO y reservoir species (i.e., N 2 O 5 , 

1075 HNO3, and CIONO 2 ) by a series of chemical processes (see 
next section) at different time scales. In general, the NO x 
deactivation is very slow in the upper stratosphere. At lower 
altitudes, however, observed HNO 3 and CIONO2 increases 
immediately after the onset of the proton forcing during the 
loso Halloween event indicate a much faster conversion. In or- 
der to assess the agreement of observed and modeled SPE- 
related odd nitrogen enhancements, we have thus compared, 
at first instance, total NO y (= NO + NO 2 + HNO3 + 2 N 2 O 5 
+ CIONO 2 + HNO 4 ) rather than NO x . Since meridional 
was redistribution is an issue (see discussion in the previous sub- 
section), we have analyzed area-weighted averages of NO y 
enhancements with respect to 26 October within 40-90°N, 
i.e., in an area that entirely covers the source region. 

As a first step, we analyze the instantaneous NO y 
1090 enhancements during the main proton forcing around 
29 October-1 November. Figure 12 shows the observed 
and modeled NO y enhancements during this period, ranging 
from a few ppbv in the middle stratosphere to several 100 
ppbv in the mesosphere. The agreement between the obser- 
1095 vations and the multi-model average, the latter providing a 
measure of the overall ability of current atmospheric models 
to reproduce SPE-related NO y increases, is reasonable, ex- 
hibiting differences below 30% in the whole altitude range. 
There is, however, a systematic overestimation (underestima- 
1100 tion) of the models around 1 hPa (above 0.3 hPa). 

The NO y underestimation of the models above 0.3 hPa 
could be related to an overestimation of NO photolysis, the 
principal NO y loss mechanism in the sunlit mesosphere. It 
has been pointed out by Minschwaner and Siskind (1993) 
1105 that absorption of solar irradiance by thermospheric NO, be-1135 
ing significantly enhanced during SPEs, has an important im- 
pact on the photolysis rates of nitric oxide in the middle at- 
mosphere. 

The systematic behavior of the NO y overestimation 
mo around 1 hPa suggests that these differences are related - abwo 
least partly - to the simulated ionization rate profile. In this 
pressure range, there are uncertainties in the modeling of 
electron precipitation at 300 keV to 5 MeV. Electrons con- 
tribute there to the total ionization within 40-90°N by ap- 
1115 proximately 15%. As the highest electron channel on POESiws 
does not provide data up to 5 MeV, the energy spectra was 
extended according to Klassen et al. (2005). However, the 
overestimated NO y at 1 hPa questions the extended spec- 
tra. In addition, the energy range of the highest electron 
1120 channel mepOeJ is not known for sure (private communica-uso 
tion, Janet Green, NOAA) and it might be smaller than the 
published 300 keV-2.5 MeV (Evans and Greer, 2000). A 
smaller energy range would give an increase to NO y pro- 
duction at 0.1 hPa. A possible overestimation of electron 
1125 ionization is also supported by the better agreement of theuss 
WACCM simulation without electrons (WACCMp) with the 



Fig. 13 . NO* production efficiency (ratio of the net NO* increase 
and the integrated initial N production) during the period of the 
main proton forcing (28 October-1 November) from box model cal- 
culations for night- and daytime conditions (solid and dashed black 
lines, respectively), assuming initial atmospheric conditions as ob- 
served by M1PAS in the polar cusp region. The following variations 
for dark conditions are also shown: a 20 K temperature increase 
(red) and decrease (blue), a factor of 2 increase of O 3 (solid green), 
and initial NO* abundances set to zero (dashed green). 

observations. Above 2 hPa, this simulation shows approx- 
imately 20% less enhanced NO y than the nominal simula- 
tion, including protons and electrons. Additional ionization 
by alpha particles, included in CAO, FinROSE, SOCOL, and 
SOCOLi contributes only by approximately 5% to the total 
ionization within 40-90°N, hence increasing the SPE-related 
NO y enhancements only marginally. 

On the other hand, there is an important spread of up to 
100% among the modeled NO y enhancements. This is sur- 
prising, given that all models use the same AIMOS ioniza- 
tion rates and, thus, rather similar modeled NO y enhance- 
ments are expected at least during the first days of the main 
proton forcing. NO y enhancements are most strongly over- 
estimated (up to 100%) by SOCOLi, SOCOL, and CAO in 
the stratosphere around 1 hPa. In the mesosphere, smallest 
NO v increases are obtained by B2dM and EMAC (up to 80% 
less than observed), while SOCOL and SOCOLi simulations 
overestimate the MIPAS enhancements by 50%. 

In order to investigate possible reasons for the spread 
among the model results, a more detailed look into the NO x 
production mechanism is required. Generally, it is assumed 
that each ion pair produces 1.25 atomic nitrogen atoms, dis- 
tributed between the electronic ground state N( 4 S) and the 
excited N( 2 D) state with a branching ratio of 0.45 and 0.55, 
respectively (Jackman et al., 2005b). The value of 1.25 
atomic nitrogen atoms per ion pair has been adapted by all 
models involved in this study, except for EMAC and SO- 
COLi. In the latter model, N production is implicitly mod- 
eled by means of the involved ion chemistry scheme. In 
EMAC, an altitude-dependent N production has been as- 
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Fig. 14. Spatial distributions of observed and modeled NO y at 1 hPa averaged of the period of the main proton forcing during 30 October- 
1 November. The average precision of MIPAS observations is also shown (upper second panel from the left). 


sumed which has been determined empirically by the adjust- 
ment of the simulations to observed NO 2 and N 2 0 abun- 
dances (Baumgaertner et al., 2010). The resulting N pro- 
1160 duction profile is slightly higher than that used by the othemeo 
models in the upper stratosphere (around 1 .5 N per ion pair) 
and considerably lower in the mesosphere (less than 0.3 N 
per ion pair) which explains to a major extent the behavior of 
the EMAC NO y enhancements compared to other models, 
lies An important source of variability in the NO y production's 
is related to the reaction paths of the produced atomic ni- 
trogen in its ground and excited states. While the reaction 
of N( 2 D) with oxygen to form NO is very fast such that 
practically all N( 2 D) is immediately converted to NO below 
iiro the thermosphere, the corresponding reaction of the nitrogermso 
ground state 

N( 4 S) + 0 2 ->N0 + O (Rl) 

is slower and highly temperature-dependent. Hence, it com-, 195 
petes with other reactions, namely: 

1175 N( 4 S) + N0->N 2 + 0 (R2) 

N( 4 S) + N0 2 ->N 2 0 + 0, (R3) 


both destroying NO x . As a consequence, only a fraction of 
the initially produced NO x remains available after the pro- 
ton forcing. This fraction depends strongly on temperature 
due to reaction Rl and to a lesser extent on the reparti- 
tioning between NO and NO 2 , driven by illumination and 
odd oxygen availability. In order to assess the sensitivity 
of the SPE-related NO x production to these parameters, we 
have integrated the relevant chemical equations for the pe- 
riod 28 October-1 November with a simple box model in- 
cluding AIMOS ionization rates and assuming initial atmo- 
spheric conditions as observed by MIPAS at 70-90°N before 
the SPE, as well as the N( 4 S)/N( 2 D) branching ratio recom- 
mended by Jackman et al. (2005b) . The modeled NO x en- 
hancements have then been compared to a similar simulation, 
but setting the rate coefficient for reactions R2 and R3 to 
zero (i.e., assuming that all initially produced NO survives). 
The ratio ofboth simulations reflects the NO x production ef- 
ficiency. It is shown in Fig. 13 for nighttime and daytime 
conditions (solid and dotted black lines, respectively), show- 
ing maximum value of 0.55-0.7 around the stratopause and 
smaller values (0.15-0.4) above and below. Reduced values 
below the stratopause are related to the background NO x : if 
initial NO x abundances are set to zero, the production effi- 
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Fig. 15. Temporal evolution of area-weighted averages of NO y changes with respect to 26 October 2003 in MIPAS observations and model 
simulations at 40-90° N, as well as differences between modeled and observed averages. The confidence of observed NO y changes (in units 
of cr) is shown in the upper second panel (from the left). Note that observations exhibiting CO abundances higher than 1 ppmv have been 
omitted in the averaging in order to exclude EEP- induced contributions to the NO y enhancements. See text for further details. 


ciency increases with pressure to values close to unity in the red and blue lines in Fig. 13). On the other hand, assuming 
lower stratosphere (see Fig. 13, dotted green line). A temper -1205 a two times higher ozone abundance results in an increase of 
ature increase (decrease) of 20 K results in an enhancement the NO x production efficiency by only a few percent, 
(reduction) of this quantity by approximately 30-50% (see 

The chemical scheme described above (including a N( 4 S) 
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Fig. 16. As Fig. 15, but showing area-weighted averages for WACCMp (proton ionization, only, left), differences between WACCMp and 
M1PAS (middle), and differences between WACCMp and WACCM (right). 


and N( 2 D) branching ratio of 0.45 and 0.55) has been em- 
ployed in most of the atmospheric models included in the 
1210 intercomparison, with some exceptions: B3dCTM and CAO 1250 
use a family approach which implies the immediate conver- 
sion of all atomic nitrogen to NO (equivalent to a ratio of 
1 in Fig. 13), explaining - at least partly - the relatively 
high NO y increases above 2 hPa in these models. Also Fin- 
1215 ROSE applies a family approach, however, in this model it isnss 
implicitly assumed that all N( 4 S) produced by ionization de- 
stroys NO via reaction R2, resulting in a net NO x production 
of 0.25 per ion pair (i.e., an altitude-independent production 
efficiency of 0.2 in Fig. 13), However, although a consid- 
1220 erably smaller NO y production is hence expected, FinROSE™ 
model results show more excess NO y than found in the ob- 
servations. EMAC uses a N( 4 S) and N( 2 D) branching ratio 
of approximately 0.2 and 0.8, respectively. Box model cal- 
culations using this atomic nitrogen branching yield weakly 
1225 altitude-dependent NO x production efficiencies of 0.6-0.8,1265 
considerably higher than the nominal efficiency of ~0.2 in 
the mesosphere. Therefore, the smaller atomic nitrogen pro- 
duction in the mesosphere applied in EMAC is partly com- 
pensated by the modified N( 4 S) and N( 2 D) branching ratio. 

1230 As shown above, temperature differences might explains 
the differences of the NO y enhancements simulated by the 
remaining models. B2dM underestimates the observed tem- 
peratures in the mesosphere by about 1 5 K, consistent with 
the relatively low NO y enhancements compared to the other 
1235 models and observations, there. In contrast, HAMMONIA275 
and SOCOLi simulations, exhibiting relatively low meso- 
spheric temperatures, show much larger NO y enhancements. 
Stratospheric temperatures are significantly overestimated by 
B2dM, CAO, HAMMONA, and SOCOLi. However, only 
1240 the latter model shows stratospheric NO y enhancements welkso 
above the model average. Thus, temperature differences 
among the models cannot be the only reason for the spread 
encountered in the modeled NO y enhancements. 

Therefore, we have looked at the spatial NO y distribu- 
1246 tion in order to investigate if the spread in the modeled NOy 1 " 85 
could also be related to dynamical effects. Figure 14 shows 
the observed and modeled NO v distributions in the upper 


stratosphere (1 hPa) averaged over the period 30 October- 
l November. The spatial extension of the modeled NO y en- 
hancements exhibits pronounced differences. In some cases, 
NO y enhancements are confined to the polar region north- 
ward of 70°N (i.e., B2dM, KASIMA, FinROSE) while in 
other cases they extend even to regions equatorwards of 50° N 
(i.e., SOCOL and SOCOLi). Taking into account that the 
spatial extension of the source region is the same in all mod- 
els, these differences must be related to transport acting on a 
very short time scale. As discussed above, SOCOLi shows 
higher NOy averages than other models with similar strato- 
spheric temperatures (i.e., HAMMON1A, B2dM). The spa- 
tial NOy distribution of SOCOLi at 1 hPa indicates strong 
wave activity resulting in a deformation of the pole-centered 
shape of the NO y distribution. Thus, it cannot be excluded 
that tropical NO y , transported into the 40-90°N region, con- 
tributed to the large NO y enhancements identified in this 
simulation. Further, the fast transport of SPE-generated NO y 
out of the source region in the SOCOLi simulations might re- 
sult in a higher net NO y production since NO x destruction 
by reactions with atomic nitrogen (R2 and R3) is then less 
efficient. 

During the following month, the SPE-induced NO y en- 
hancements were transported downwards with the merid- 
ional circulation, forming a NO y layer around 45 km at the 
end of November (Lopez-Puertas et al., 2005a). At the same 
time, NO x , generated by continuous EEP in the lower ther- 
mosphere, reached the upper stratosphere and began to merge 
with the upper part of the SPE-induced layer (see also Fig. 1 ). 
This behavior is not reproduced by the models since low and 
mid-energy EEP is not included in the majority of the mod- 
els. In order to facilitate the comparison of observed and 
modeled SPE-induced NO y enhancements in the following 
month after the proton forcing, we have excluded those parts 
of all observed and modeled NO y profiles where M1PAS CO 
abundances were higher than 1 ppmv. This value has been 
chosen such that the major fraction of EEP-induced NO y en- 
hancements has been filtered out without removing too many 
MIPAS locations, particularly at higher altitudes. 

Figure 1 5 shows the temporal evolution of the observed 
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Fig. 17. Spatial distributions of observed and modeled NO y at 2 hPa averaged of the period 20-27 November. 


and modeled NO y enhancements (related to the SPE, only) 
with respect to 26 October within 40-90°N for the following 
i29o month. While the magnitude of the enhancements is gen- 
erally larger than in the observations and further shows a 
significant spread related to the differences in the NO y pro- 325 
duction during the proton forcing (see discussion above), the 
observed evolution of the SPE-induced NO x layer is well re- 
1295 produced by all models in terms of vertical distribution and 
relative vmr decrease. A more detailed look into the tempo- 
ral NO y evolution of individual models shows that smaller^ 
fluctuations can be attributed to dynamical variability. 

It is interesting to notice that the WACOM simulation 
1300 without electron-induced ionization yields a much better 
agreement with the observations than the nominal simulation 
throughout the period under investigation (see Figure 16).i335 
Additional NO y buildup related to electron-induced ioniza- 
tion is even more pronounced during the second event (4- 
1305 5 November) below 0.4 hPa compared to the main proton 

forcing (see right panel of this Figure). This excess produc- 
tion during the second event dominates the NO y overesti- 1340 
mation of 5- 10 ppbv encountered in the nominal WACCM 
simulation during the following weeks. 

1310 The meridional distributions of the observed and mod- 
eled NO y enhancements exhibit important differences to- 
wards the end of November (see Figure 17, showing NO y i345 
distributions at 2 hPa averaged over 20-27 November). The 
observed and modeled latitudinal gradients correlate well 
1315 with the meridional CFI4 anomalies (see Fig. 9), highlight- 
ing the important role of mixing and large-scale transport. 
The meridional redistribution of the SPE-induced NO y en -1350 
hancements, particularly the transport out of the polar night 
region, has important implications on the NO y repartitioning 
1320 which is to a major part driven by photochemistry (see next 


subsection). 

An interesting detail of the observed evolution of SPE- 
induced NO y enhancements (Fig. 15, upper left panel) is 
the appearance of several ’’spikes” at mesospheric altitudes, 
which are temporally correlated with peaks in the ionization 
related to high energy (>300 keV) electron precipitation, the 
most pronounced event occurring on 21 November. Ion- 
ization by high-energetic electrons is included in the mod- 
els which, however, do not reproduce such sudden NO y in- 
creases. It is therefore unlikely the observed mesospheric 
NO y peaks are related to in situ production by EEP associ- 
ated to the Halloween event. Instead, they could be related 
to residual contributions of descending NO x from the up- 
per mesosphere which have not completely been filtered out. 
It should be noted that the observed CO temporal evolution 
(see Figure 1 1 ) indicates particularly strong descent around 
20 November. 

A fraction of the NO x deactivation by reaction with atomic 
nitrogen during the proton forcing discussed above occurred 
via reaction R3, giving rise for the buildup of N 2 O. Upper 
stratospheric and mesospheric nitrous oxide increases up to 
7 ppbv have been observed by MIPAS during the Halloween 
SPE and have been attributed to this reaction channel (Funke 
et al., 2008). Reasonable agreement with CMAM model cal- 
culations has been obtained by assuming that only half of the 
products of reaction R3 is N 2 0 and O, while the other half is 
N 2 and 0 2 . Figure 1 8 shows the observed and modeled N 2 0 
zonal mean enhancements averaged over the period of the 
main proton forcing (29 -31 October). Except for FinROSE 
and B2dM, which do not include the reaction channel R3, 
N 2 0 increases are simulated by all models. The observed en- 
hancements, however, are generally overestimated by a fac- 
tor 2 to 10, except for EMAC which shows smaller N 2 0 in- 
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Fig. 18. Zonal mean NjO changes with respect to 26 October 2003 in MIPAS observations and model simulations averaged over the period 
29-31 October. The significance of observed NeO changes (in units of a) is shown in the upper second panel (from the left). 


creases than observed by MIPAS. In the latter model, this , 375 
can be clearly attributed to the modified N( 4 S) and N( 2 D) 
branching ratio (see discussion above). Except for WACCM, 
the remaining models do not include the additional reaction 
channel of R3, responsible for the formation of N 2 and O 2 , 
which has been included in the CMAM simulations (Funke 
et al., 2008). But even when taking into account a reduction 
by a factor of 2 of the simulated enhancements, these mod-, 36 o 
els tend to overestimate the observations and further show 
a significant spread among the individual results. As in the 
case of NO y , also the total SPE-induced N 2 O production 
depends on temperature, NO x partitioning, and dynamical 
redistribution. However, a dominant relationship of none of 
these quantities with the differences of the magnitudes of the 1385 
modeled N 2 0 increases can be established. 

6.3 Repartitioning of nitrogen species 

1390 

After having assessed the observed and modeled total NO y 
and N 2 0 enhancements generated by the Halloween event, 
we analyze in this subsection the repartitioning of initially 
produced nitric oxide into other NO y species in the aftennath 
ofthe SPE. 1395 


6.3.1 NO x 

The conversion of the excess NO generated by the proton 
forcing into NO 2 acts on a very short timescale (seconds to 
minutes) and is controlled at dark conditions by the reactions 

NO + O 3 — t-N0 2 ~l-0 2 (R4) 

N0 2 +0-s-N0 + 0 2 , (R5) 

giving rise to a N0 2 /NO x ratio close to one in the strato- 
sphere, but decreasing in the mesosphere due to the availabil- 
ity of atomic oxygen. This decrease occurs at higher altitudes 
in the polar night region compared to midlatitudes. Figure 19 
shows the observed and modeled nighttime N0 2 / NO x ratios 
averaged over the initial SPE period. The observed decrease 
of this ratio above 0.3 hPa at midlatitudes and 0.1 hPa in 
the polar region is generally well reproduced by the models 
which resolve the mesosphere, except for B2dM and EM AC, 
which both overestimate the polar NO 2 fraction at these al- 
titudes. The higher mesospheric N0 2 abundances in these 
two models might be related to lower atomic oxygen con- 
centrations at high altitudes and/or less efficient mixing be- 
tween polar night and illuminated regions. At sunlit condi- 
tions, photolysis of N0 2 and higher atomic oxygen abun- 
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Fig. 19. Zonal mean nighttime NCh/NOx ratios averaged over the period 28 October-15 November in M1PAS observations and model 
simulations. 


dances shift the N0 2 /N0 x ratio to lower values compared to 
dark conditions. Figure 20 shows the observed and mod- 
eled daytime ratios. The observed values are well repro- 
duced by B2dM, SOCOL, and SOCOLi, while other models 
moo tend to overestimate the polar upper stratospheric and meso -1420 
spheric NO 2 fraction close to the terminator. These differ- 
ences in the NO x partitioning among the models and obser- 
vations highlight the difficulties in drawing conclusions on 
the SPE-induced total NO x enhancements from the compar- 
Mos ison if only one of its components is considered. 

1425 

6.3.2 N 2 O 5 and HNO 3 

In the stratosphere, excess NO x produced by the proton forc- 
ing, is slowly converted into N 2 Os in the following weeks 
after the SPE by the reactions 1430 

i 4 io NO 2 + O 3 — 4 NO 3 |- O 2 1R6) 

N0 2 + N0 3 + M-»N 2 0 5 + M. (R7) 

The rate-limiting reaction of this conversion is R 6 whichios 
exhibits a strong temperature dependence. N 2 O 5 enhance- 
ments around 1-2 ppbv, appearing 10-1 5 days after the Hal- 
1415 loween event, have been observed by MIPAS around 70- 


90°N and have been attributed to the repartitioning of SPE- 
induced excess NO x (Lopez-Puertas et al., 2005b). This con- 
version is further accelerated in the course of November by 
the growth of the polar night region, reducing the efficiency 
of N 2 O 5 losses by photolysis. 

Figure 21 shows the observed and modeled temporal evo- 
lutions of the N 2 Os enhancements with respect to 26 Octo- 
ber averaged over 70-90°N. A stratospheric N 2 O 5 buildup, 
being most pronounced in the second half of November, is 
simulated by all models, qualitatively reproducing the ob- 
served behavior. The modeled N 2 O 5 increases are, however, 
generally overestimated (except for KASIMA) and exhibit a 
wide spread among the models. Taking into account that the 
magnitude of the N 2 O 5 increase depends on various factors 
such as NO x availability, temperature, ozone abundances, 
and the efficiency of N 2 O 5 -+ HNO 3 conversion (see be- 
low), a large spread of the model results is expected. B2dM 
and EMAC, however, overestimate the observed N 2 Or> in- 
creases by factors of 4 and 6 , respectively. While in the case 
of B2dM the extraordinarily high N 2 Og amounts can be ex- 
plained by the very pole-centered distribution of the precur- 
sor NO x , implying insignificant photochemical losses in the 
source region (see also discussion below), the reason for the 
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Fig. 20. Zonal mean daytime N02/NO x ratios averaged over the period 28 October-15 November in MIPAS observations and model 
simulations. 


unreasonably high N 2 O 5 abundances of tip to 12 ppbv in the 
mo case of EMAC is still under investigation. N 2 O 5 enhance- 
ments simulated by CAO until 4 November are likely to be 
caused by seasonal variations rather than by the SPE. «5 

Two distinct HNO 3 enhancements were observed by MI- 
PAS in the aftermath of the Halloween SPE (Lopez-Puertas 
1445 et al., 2005b). The first one, reaching vmrs around 2 ppbv, 
occurred immediately after the SPEs at altiudes above 40 km 
and has been initially attributed to the gas-phase reaction 
NO 2 + OH + M -t HNO 3 + M. Verronen et al. (2008), how- 
ever, have shown that the instantaneous HNO 3 increase af- 
1450 ter the proton forcing can only be reproduced by model cal- 
culations including ion-ion recombination between NO, 3 ~ 
and H + cluster ions. The second enhancement of 1 -5 ppbv 

1475 

started around 10 November and lasted until the end of De- 
cember. Also in this case, attempts to reproduce the magni- 
1455 tude of the observed increases by model calculations includ- 
ing gas phase chemistry only, have failed (Jackman et al., 
2008). 

1480 

Figure 22 shows the observed and modeled temporal evo- 
lutions of the HNO 3 enhancements with respect to 26 Octo- 
i46o ber averaged over 70-90°N. Consistent with previous find- 
ings, the first instantaneous enhancement is considerably un- 


derestimated by all models, except FinROSE, which includes 
the ion chemistry proposed by Verronen et al. (2008). This 
model, however, overestimates the observed increases by up 
to a factor of 3. The overestimation below. 50 km is surpris- 
ing, given that 1-D simulations with the Sodankyla Ion and 
Neutral Chemistry (SIC) model which includes the ion-ion 
recombination were found to be in good agreement with the 
same MIPAS dataset at these altitudes. It should be noted, 
however, that different ionization rates have been used in 
the SIC calculations of Verronen et al. (2008), and -probably 
even more important - that FinROSE uses a parameterization 
of the ion-ion recombination included in the full ion chem- 
istry scheme of the SIC model. 

The second enhancement, occurring around 1 5 November 
at l-2hPa, is only reproduced by KASIMA, however, over- 
estimating the observed increases by a factor of 3. Contrary 
to other models, KASIMA simulations account for HNO 3 
formation via water cluster ions (Bohringer et al, 1983) 
combined with heterogeneous reactions on sulfate aerosols 
by means of a parameterization provided by de Zafra and 
Smyshlyaev (2001 ). At lower altitudes (i.e., below lOhPa), 
midterm HNO 3 increases are visible in the observations, 
as well as in the B2dM, B3dCTM, EMAC, FinROSE, and 
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Fig. 21. Temporal evolution of area-weighted averages of NoOr, changes with respect to 26 October 2003 in M1PAS observations and 
model simulations at 70-90°N. as well as differences between modeled and observed averages. Solid contour lines reflect 1 ppbv steps. The 
confidence of observed N 2 O 5 changes (in units of u) is shown in the upper second panel (from the left). 


mss WACCM model results. These increases are not related to 
the SPE and can be explained by seasonal variations. 

In order to assess the repartitioning of the main NOyisio 
species towards the end of November in a more quantita- 
tively way, we have analyzed their relative contributions to 
mm the total NO y . This is necessary because of the encoun- 
tered differences in the total amount of SPE-induced excess 
NO y among the different models and the observations. Duei 5is 
to the observed conversion of N 2 O 5 into HNO 3 , we have 
looked, as a first step, at the relative contribution of the sum 
1495 of both reservoir species to NOy. Observed and simulated 
zonal mean ( 2 N 2 05 +F 1 N 03 )/N 0 y ratios, averaged over the 
period 15-30 November, are shown in Fig. 23. The observed®,) 
ratio of 0.28 at the peak height (~0.2hPa, see Fig. 15) of 
the NOy enhancements in late November (indicated by a 
1500 black line in Fig. 23) is very well reproduced by all mod- 
els, except B2dM and EMAC. As discussed above, the dis- 
agreement found in these models is produced by too effi- 
cient buildup of N 2 Or> (see Fig. 21). The differences ofezs 
the ( 2 N 205 +HN 0 s)/N 0 y ratio in the observations and the 
1505 B2dM simulations are, however, much less pronounced than 
those encountered in the absolute N 2 0 5 abundances: while 


B2dM N 2 O 5 exceeds the observed amounts by a factor of 
4, the modeled 2 N 2 O 5 +HNO 3 contribution to NO y at its 
peak height is around 40%, exceeding the observed contri- 
bution by only a factor of 0.5. The N 2 O 5 overestimation 
in this model is hence mainly related to the higher amounts 
and more pole-centered distribution of the precursor NO 2 . In 
contrast to B2dM , EMAC simulations obtain more than 90% 
of the available NO y at its peak altitude in the formofN 2 05 . 
This contribution decreases with altitude, but still exceeds 
30% in the mesosphere. Other models show, in some cases, a 
minor overestimation of the reservoir species fraction which 
can be partly explained by differences in the modeled tem- 
peratures and ozone abundances, controlling the efficiency 
of reaction R 6 . 

The repartitioning between HNO 3 and N 2 O 5 has been as- 
sessed by comparing the observed and modeled zonal mean 
HN 03 /( 2 N 2 0 5 +HN 03 ) ratios averaged over the period 1 5- 
30 November (Fig. 24). As expected, the observed ratio 
is strongly underestimated above approximately lOhPa by- 
all models, except KASIMA. The qualtitative agreement of 
KASIMA simulations and MIPAS observations is very good, 
particularly regarding the vertical shape of this ratio. The 
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Fig. 22. Temporal evolution of area-weighted averages of HN0 3 changes with respect to 26 October 2003 in MIPAS observations and 
model simulations at 70-90°N, as well as differences between modeled and observed averages. Solid contour lines reflect 1 ppbv steps. The 
significance of observed HNO 3 changes (in units of a) is shown in the upper second panel (from the left). 
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Fig. 23. Zonal mean ( 2 N 20 S +HN 03 )/N 0 y ratios averaged over the period 1 5-30 November in MIPAS observations and model simulations. 
The MIPAS ratio of 0.28, encountered at altitude of the maximum of the SPE-induced NO v layer, is indicated by a black line. 
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Fig. 24. Zonal mean HNO.'j/( 2 N 2 05 +HNC> 3 ) ratios averaged over the period 15-30 November in M1PAS observations and model simula- 
tions. Regions with observed vmrs of 2 N 2 O 5 +HNO 3 smaller than 0. i ppbv have been omitted. 


modeled ratio, however, exhibits a positive bias of 0.2 with 
1530 respect to the observations, most pronounced in the polar re-isec 
gion. We conclude that the HNO 3 formation via water clus- 
ter ions and/or heterogeneous reactions on sulfate aerosols, 
both included in KASIMA by means of the parameterization 
of de Zafra and Smyshlyaev (2001), is the responsible mech- 
1535 anism for the observed HN0 3 enhancements in late Novem - !565 
ber. However, some further work is required to adjust the 
parameterization quantitatively to the measurements. 

6.3.3 Minor NO y species 

Also minor NO y species were found to be enhanced in the af- 1 5?0 
1540 termath of the Halloween SPE due to the repartitioning of ini- 
tially produced NO x . Lopez-Puertas et al. (2005b) reported 
CIONO 2 enhancements up to 0.4 ppbv a few days after the 
proton forcing from MIPAS observations. These observa- 
tions are compared to the model simulations in Section 6.6, 1575 
1545 together with observations of other chlorine species. MI- 
PAS has also observed enhanced HNO 4 during the first days 
of the Halloween SPE which have not been reported so far. 
These increases can be attributed to the termolecular reaction 

1580 

H0 2 + N0 2 + M->HN0 4 + M. (R 8 ) 

1550 Since at the altitude of the HNO 4 enhancements (around 2- 
3 hPa) SPE-related increases of the precursor NO 2 are rel- 
atively small (~ 2 ppbv) compared to the background NO 2 
abundance, the observed HNO 4 changes are mainly driven 1585 
by enhanced HO 2 abundances, and hence, represent an indi- 
1555 cator of SPE-generated HO x in the middle stratosphere. At 
dark conditions. HO x is in steady state even during a SPE, 
and its abundance is hence directly proportional to atmo- 
spheric ionization. Stratospheric HNO 4 is destroyed duringsao 


the day by photolysis and by reaction with OH. Nighttime 
losses are negligible under quiescent conditions, and even 
during SPEs, OH-related HNO 4 destruction is small com- 
pared to its production via reaction R 8 . 

Due to problems with the gain calibration, particularly 
affecting this species (see discussion in Section 2), we re- 
strict our analysis to data from the gain calibration period 
28 October - 5 November. Figure 25 shows the observed 
and modeled zonal mean distributions of HNO 4 vmrs dur- 
ing the first four days of the proton forcing (29 October- 
1 November). Model results for pre-SPE conditions (26 Oc- 
tober) are also shown. HNO 4 model output is not available 
from CAO, EMAC, HAMMONIA and FinROSE. Polar up- 
per stratospheric enhancements of up to 0. 1 8 ppbv are vis- 
ible in the observations on 29 October, decreasing until 1 
November by about 20%. The HNO 4 enhancements are also 
simulated by the models in the first days of the SPE, how- 
ever, generally overestimating the measurements. The over- 
estimation is most pronounced in the B2dM, B3dCTM, and 
WACOM simulations (a factor 2-3), while SOCOLi shows 
smaller HNO 4 increases. The HNO 4 peak height is located 
at somewhat lower altitudes in SOCOLi which might be re- 
lated to the relatively high abundances inside the ambient 
HNO 4 layer around 5 hPa. Both, the SPE-related and ambi- 
ent peaks can not be vertically resolved and merge together 
after the application of MIPAS averaging kernels. The mod- 
erate decrease of HNO 4 in the following days is qualitatively 
reproduced by all models except B2dM. In this particular 
model, the HNO 4 enhancements are confined to the polar 
night region, hence experiencing less photochemical losses. 

The differences in the magnitude of the HNO 4 enhance- 
ments in the simulations and the observations can partially 
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Fig. 25. Observed and modeled zonal mean HNO 4 vmrs for pre-SPE conditions (26 October) and during the main proton forcing 
(29 October-- 1 November). Solid contour lines reflect 0.1 ppbv steps. Note that MIPAS observations from 26 October have been omit- 
ted due to gain calibration problems. 


be explained by differences in the abundances of the precur- 
sor N0 2 . During the main proton forcing, modeled NO 2 
abundances at 70-90 c N at the HNO 4 peak height are on av- 
erage 50-100% higher than the observed ones (not shown). 

Isas The NO 2 overestimation in the models is mainly related to eon 


a weaker degree of denoxification compared to the observa- 
tions already before the SPE event. Aditionally, differences 
in the HO 2 availability might also play an important role 
in explaining the behavior of modeled HNO 4 At the peak 
height of the ITNO 4 enhancements and in the absence of sun- 
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Fig. 26. Temporal evolution of area-weighted averages of relative O 3 changes with respect to 26 October 2003 in MIPAS observations and 
model simulations at 70-90°N, as well as differences between modeled and observed averages. The significance of observed O 3 changes (in 
units of a) is shown in the upper second panel (from the left). 


light, HO2 is the dominant HO x compound in the presence 
of proton forcing. Since FIO x productions by SPE-induced 
ionization is not expected to differ significantly among the 
models, different HOv abundances are most likely related to 
differences in the EIO x losses. These are dominated by the 


cannibalistic reaction 


In this sense, the HO x partitioning plays a crucial role in the 
HO2 availability at these altitudes despite the approximately 
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Fig. 27. Area-weighted averages (70-90°N) of observed and mod- 
eled relative O 3 changes during 29 October— 4 November (left) with 
respect to the abundances on 26 October (right). Thick solid and 
dashed lines represent multi-model mean average and MIPAS ob- 
servations, respectively. WACCMp denotes the WACOM simula- 
tion including proton ionization, only (excluded from the multi- 
model mean). reso 

1610 20 times smaller OH concentrations compared to HO 2 . 

It might also be possible that the HO x partitioning is af- 
fected by ion chemistry. Several ion chemistry reactions are™ 
known which transfer H into OH~, and therefore might act 
as a sink of HO 2 ; one reaction is known which transfers OH 
1615 into HO 2 . While it is beyond the scope of this investigation 
to determine whether these reactions really significantly af- 
fect the partitioning between odd hydrogen species, it mighte7o 
be worthwhile to investigate this point in the future. 

6.4 Ozone loss 

1620 One of the most important aspects of the model-data inter -1675 
comparison of SPE-induced composition changes during the 
Halloween event is the evaluation of the ability of the mod- 
els to reproduce the observed ozone destruction caused by 
acceleration of catalytic HO x and NO x cycles. SPE-induced 
1625 ozone losses have been observed by a variety of space-bome 6 8o 
instruments during several of the stronger events of the past 
two solar cycles (see Jackman et al., 2000, for a review). 
Two different types of ozone destruction could be distin- 
guished: HO x -related short-time losses, acting principally in 
i63o the mesosphere during the event itself, and NO x -related mid-™ 
term losses in the stratosphere which can last up to several 
months in the polar winter atmosphere. Such a behavior has 
also been observed by MIPAS in the aftermath of the Hal- 
loween event. Lopez-Puertas et al. (2005a) reported HO x - 
1635 driven mesospheric ozone losses up to 70% and NOx-drivem® 
stratospheric losses of around 30%, the latter lasting for more 
than 2 weeks in the Northern hemisphere. 

Figure 26 shows the observed and modeled temporal evo- 
lutions of the relative O 3 changes with respect to 26 Octo- 
1640 her, averaged over 70-90°N. The mesospheric ozone losses™ 


above 0.3 hPa, which exhibit two distinct peaks related to the 
main proton events on 29 October and 4 November, are well 
reproduced by most of the models. Also the stratospheric O 3 
losses during the following month, peaking around 1 hPa, are 
qualitatively reproduced by the simulations, however, with a 
more pronounced spread of the model results. This is not 
surprising since these losses are driven by NO x which ex- 
hibits important differences between the models, particularly 
during the second half of November (see sections 6.2 and 
6.3). Further, NO x -induced ozone loss is driven by reac- 
tion R4 which is very sensitive to temperature differences. 
The midterm evolution in the mesosphere is characterized by 
ozone buildup which is related to seasonal variations (sum- 
mer to winter transition), and which is generally more pro- 
nounced in the model simulations compared to the observa- 
tions. 

In order to assess observed and modeled short-term ozone 
depletion in a more quantitative way, we have compared pro- 
files of relative ozone changes at 70-90°N, averaged over 
the period of the main proton forcing (28 October-4 Novem- 
ber), in Figure 27. The agreement between observations and 
the multi-model average is excellent in the mesosphere, in- 
dicating a very good overall ability of the models to repro- 
duce HO x -related ozone losses under SPE conditions. Also, 
the models themselves agree reasonably well in this altitude 
range, except for B2dM. In the stratosphere, where NOx- 
related losses are dominant, the agreement between the mod- 
els is worse, though the model average is very close to the ob- 
servations within 5%. Ozone depletion around 1 hPa is over- 
estimated by EMAC and B3dCTM. CAO and SOCOL re- 
sults indicate a somewhat smaller ozone loss throughout the 
stratosphere. WACCM simulations performed with and with- 
out electron-induced ionization (WACCM and WACCMp, 
respectively in Figure 27) indicate an additional ozone loss 
induced by electrons in the order of 5% above 2 hPa. It is 
interesting to notice that the agreement of the WACCM sim- 
ulations with the observation is better when excluding the 
electron contribution. 

Figure 28 shows the corresponding zonal mean distribu- 
tions. Observed mesospheric losses extend to around 60 C N 
in consonance with the expected cut-off latitude of proton 
precipitation. This latitudinal distribution is well reproduced 
by the models. B2dM shows a mesospheric ozone buildup 
poleward of 80°N related to seasonal changes, which over- 
compensates IIO x -related losses at these particular latitudes. 
This behavior, which can be attributed to deficient meridional 
mixing in the polar region, give rise for the apparent underes- 
timation of mesospheric ozone losses of B2dM in Figure 27. 
Oscillations encountered in the CAO ozone changes above 
1 hPa at 40-50^ are related to the background O 3 and are 
not caused by the SPE. 

The latitudinal extension of observed and modeled strato- 
spheric ozone losses around 1 hPa correlates well with the 
area of NO x increases shown in Fig. 14. In this altitude re- 
gion, ozone depletion is restricted to latitudes poleward of 
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Fig. 28. Zonal mean relative Os changes with respect to 26 October averaged over the period 28 October-4 November in M1PAS observations 
and model simulations. Solid contour lines reflect 20% steps. The significance of the observations is also shown (second top panel from the 
left, in units of a). 



1710 


Fig. 29. Area-weighted averages (70-90°N) of observed and mod- 
eled relative O 3 changes with respect to 26 October during 16- 
26 November . Thick solid and dashed lines represent model mean 
average and MIPAS observations, respectively, WACCMp denotes 
the WACOM simulation including proton ionization, only (excluded 1715 
from the multi-model mean). 


70°N. It is interesting to notice that B2dM simulations show 
no NO x -induced ozone loss in the upper stratospheric polan? 2 o 


night region, in contrast to the observations and other mod- 
els, Indeed, the NO x catalytic cycle is expected to be ineffi- 
cient at dark conditions since NO 2 is not reconverted to NO. 
Strong mixing is hence required in order to obtain a homo- 
geneous ozone distribution in the polar stratosphere as found 
in the observations. Ozone increases occur in the SOCOL 
simulation below 1 hPa which can be related to intrusions of 
mid-latitude air into the polar region, over-compensating the 
SPE-induced ozone losses. 

Figure 29 shows profiles of stratospheric mid-term 
ozone changes at 70-90°N, averaged over the period 16- 
26 November. As expected, modeled ozone depletions have 
a larger spread than during the main proton forcing, rang- 
ing from 10 to 50% at the peak height. The model average, 
however, is in very good agreement with the observed deple- 
tion of 30% at 1-2 hPa. Only minor differences of 5% are 
found at its maximum. WACOM simulations with and with- 
out electron-induced ionization, however, suggest that these 
remaining differences could be significantly reduced when 
excluding the electron contribution to atmospheric ioniza- 
tion. 

Figure 30 shows the corresponding zonal mean distribu- 
tions. Generally, the magnitude of the stratospheric ozone 
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Fig. 30. Zonal mean relative Oa changes with respect to 26 October averaged over the period 16-26 November in MIPAS observations and 
model simulations. Solid contour lines reflect 20% steps. 


loss at 70-90°N is anti-correlated to its latitudinal extension 
which, in turn, is linked to the spatial distribution of the SPE- 
related NO x layer (see Fig. 17). Meridional redistribution 
is hence a key factor for explaining the differences in the 
t725 modeled ozone depletions shown in Fig. 29. In particular; 755 
SOCOL simulations indicate strong meridional distribution 
around 1 hPa, resulting in polar higher ozone abundances 
than in the other models, despite of the relatively high NO y 
availability shown in Fig. 17. There, NO x -driven ozone loss 
1730 is partly compensated by in-mixing of Oa-rich air-masses 
from lower latitudes. 17eo 

Observed mesospheric ozone changes in late November 
are characterized by a pronounced increase around the polar 
night terminator which is related to the buildup of the third 
1735 ozone maximum (Marsh et ah, 2001). This rapid buildup is 
responsible for the short lifetime of HO x -related ozone de- 
pletion at these altitudes. Only in the polar night region, re- 765 
duced ozone abundances are found until the end of Novem- 
ber. This behavior is well reproduced by EMAC, KAS1MA, 

1740 WACOM and, to a lesser extent, HAMMONIA. B2dM be- 
haves in an opposite way. 

In summary, SPE-related short- and midterm ozone™ 
changes are well reproduced by the atmospheric models on 
average, though individual model results can vary signifi- 
es cantly due to differences in dynamical and meteorological 
background conditions. The good agreement between mod- 
els and observations in the mesosphere can be interpreted775 
as a verification of the parameterization of HO x production 
by atmospheric ionization included in the models. On the 
1750 other hand, WACCM simulations with and without electron- 
induced ionization suggest that the agreement between mod- 
els and observations could be even improved when excludingTso 


the electron contribution to atmospheric ionization. 

6.5 Enhancements of H 2 O 2 

MIPAS observed H 2 O 2 increases of short duration imme- 
diately after the Halloween SPE in polar night stratosphere. 
H 2 O 2 is formed by the reaction 

H0 2 + HO 2 — t H 2 O 2 + O 2 (R10) 

and is hence - together with HNO4 - an indicator for SPE- 
generated HO x in the stratosphere. During daytime, it is 
photolyzed within several hours to a day, or destroyed by the 
reaction 

H 2 0 2 + 0H->H 2 0+H0 2 . (Rll) 

Chemical nighttime losses are negligible at quiescent con- 
ditions. The availability of OH during periods of proton forc- 
ing allows for H 2 0 2 destruction also at night. These losses, 
however, are most important above the stratopause. In the 
dark stratosphere, reaction R1 1 is expected to deplete H 2 0 2 
by less than 10%. Therefore, observed H 2 0 2 increases are 
primarily driven by the production mechanism R10. 

Model output of H 2 0 2 is available from B2dM, B3dCTM, 
FinROSE, HAMMONIA, KASIMA, and WACCM. Fig- 
ure 31 shows observed and modeled zonal mean H 2 0 2 
changes during the period of the main SPEs (28 October- 
4 November). The observed increases of up to 0. 1 ppbv are 
considerably overestimated by the simulations by a factor 
of 4-7. This huge difference between observed and mod- 
eled H 2 0 2 increases can hardly be explained by a possible 
overestimation of the ionization rates by a factor of 1 .2-2, as 
suggested from the comparison of NO y increases. Although 
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Fig. 31. Zonal mean H 2 0 2 changes with respect to 26 October averaged over the period 28 October-4 November in M1PAS observations 
and model simulations. Solid contour lines reflect 0. 1 ppbv steps. The significance of observed H 2 0 2 changes (in units of tr) is shown in the 
upper second panel (from the left). 


H 2 O 2 production depends quadratically on H0 2 , total HO x i8io 
scales with the square root of the ionization rate due to re- 
action R9, being the principal chemical loss mechanism at 
nighttime. Thus, 4—7 times lower ionization rates would be 
required in order to reduce modeled H 2 O 2 increases to the 
observed values. As already mentioned in the discussion of 
HNO 4 enhancements, the availability of HO 2 during nightieis 
time SPE conditions is largely controlled by the HO x par- 
titioning. At the peak height of the H 2 O 2 increases (0.5— 

1 hPa), this dependence is even more pronounced than at the 
pressure levels of the HNO 4 enhancements (2-3 hPa) due to 
the increasing OH contribution to HO x with altitude. Thus,i 8 2 o 
the disagreement of observed and simulated H 2 O 2 hints at 
an underestimation of the OH/FIO 2 ratio in the upper polar 
stratosphere during the proton forcing. Alternatively, H 2 O 2 
formation by reaction RIO might be significantly overesti- 
mated in the models. 1825 


Meridional transport to illuminated latitudes, where H 2 O 2 
is photochemically destroyed, could also affect the magni- 
tude of the SPE-related enhancements. H 2 O 2 distributions 
simulated by B2dM, which has a very strong mixing bar- 
rier, might hence experience less photochemical losses than™ 
in other models. In fact, B2dM enhancements are more con- 
fined to the polar night region. Other models, however, show 
a similar meridional distribution as observed by MIPAS. It is 
thus unlikely, that differences in the efficiency of photochem- 
ical losses related to transport can explain the pronounced 8 ® 
differences between observed and modeled H 2 O 2 enhance- 
ments. 


6.6 Enhancements and repartitioning of chlorine 
species 

Enhancements of the chlorine species CIO, HOC1, and 
CIONO 2 have been detected by MIPAS in the aftermath of 
the Halloween SPE in the NH polar stratosphere (von Clar- 
mann et al., 2005; Lopez-Puertas et al., 2005b). Short-term 
CIO and HOC1 increases of the order of 0.2 ppbv occurred 
immediately after the onset of the proton forcing on 29 Oc- 
tober. CIONO 2 increases up to 0.4 ppbv appeared approx- 
imately 2 days later, remaining in the stratosphere for sev- 
eral weeeks. SPE-related HOC1 increases have also been 
observed by MLS on Aura during the January 2005 proton 
event (Damiani et al., 2008). These enhancements were ac- 
companied by a HC1 decrease of similar magnitude, thus 
clearly demonstrating SPE-induced chlorine activation. The 
conversion of HC1 to active species occurred in presence of 
enhanced OH via the reaction 

HC1+0H->C1+H 2 0. (R12) 

HC1 can also be incorporated into negative ions, from which 
chlorine is released mainly in the form of atomic chlorine or 
chlorine monoxide. There are also reverse reactions releas- 
ing HC1, however, it has been shown in a recent publication 
(Winkler et al., 2009) that during large solar proton events, 
chlorine activation dominates, and negative ion reactions can 
act as a significant sink of HC1, and a source of active chlo- 
rine. Atomic chlorine is rapidly converted to CIO by 

Cl + 0 :i -> CTO + 0 2 . (R13) 
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Fig. 32. Zonal mean CIO changes during the main proton forcing (29-3 1 October) with respect to 26-27 October in MIPAS observations 
and model simulations. Solid contour lines reflect 0.05 ppbv steps. The significance of observed CIO changes (in units of a) is shown in the 
upper second panel (from the left). 


In the polar night stratosphere, where SPE-generated HO x is 
dominated by HO 2 , CIO is further converted to HOC1: 

C10 + H0 2 -3H0Cl + 0 2 . (R14f° 

wo The chemical lifetime of nighttime HOC1 is very long below 
the stratopause. Above and at sunlit conditions, HOC1 is 
removed by the reaction 

HOCI+ OH -f C10+H 2 0 (R15)' 865 

and rapid photo-dissociation during daytime. Due to the high 
1845 H0 2 /0H ratio in the nighttime stratosphere under SPE con- 
ditions, active chlorine is expected to be mainly in the form 
of HOC1 due to reaction R14. Hence, ambient CIO should,,,,,, 
be reduced in contrast to the results obtained by von Clar- 
mann et al. (2005) from MIPAS observations during the Hal- 
i85o loween SPE. On the other hand, OH is the dominant HO x 
constituent during daytime and HOC1 is quickly photolyzed 
even at high solar zenith angles. CIO enhancements might e75 
hence occur in the illuminated stratosphere, if SPE-related 
HO x increases were well above the background concentra- 
1855 tion. In fact, the CIO enhancements observed by MIPAS on 
29-30 October 2003 (von Clarmann et al., 2005) took place 
outside the polar night region. 


The mid-term evolution of polar ambient CIO during the 
period of the Halloween event is characterized by a continu- 
ous decrease related to seasonal variations (see Fig. 3) which 
makes the analysis of SPE-induced changes on a longer 
timescale difficult. Therefore, we restrict our analysis to the 
period of the main proton event on 29-3 1 October. Figure 32 
shows observed and modeled changes of the CTO zonal mean 
distribution, averaged over these days, with respect to 26- 
27 October. CTO increases of ~0.1ppbv are found in the 
MIPAS observations at latitudes around 60°N in qualitative 
agreement with the previous analysis of von Clarmann et al. 
(2005). These enhancements are reproduced by none of the 
models. Evidently, simulated daytime HO x increases are 
too small compared to the ambient HO x abundances to alter 
noticeably the CIO availability. Although the observed en- 
hancements are significant at the 2er-level with respect to the 
average measurement precision (see Fig. 32, second panel), 
this important difference between the observations and the 
simulations should be carefully interpreted due to a possi- 
ble systematic bias related to gain calibration errors in the 
measurements (see Section 2), particularly because the ob- 
served CTO change has been calculated from temporal aver- 
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Fig. 33. Zonal mean CIO vmrs before the SPE 26-27 October n MiPAS observations and model simulations. Solid contour lines reflect 
0. 1 ppbv steps. 


i 860 ages belonging to different gain calibration periods. In the 
polar night region, both, observations and models show a 
CIO decrease. The observed CIO reduction of up to 0.2 ppbvisos 
is considerably underestimated by the simulations, except for 
CAO. The latter model overestimates the CIO reduction by 
lass approximately a factor of 10. The unreasonably large CIO 
depletion in CAO is related to a high CIO availability be- 
fore the SPE and goes along with a C10N0 2 buildup of a 
similar magnitude (see below). The reason for the higher 19 ’ 0 
background clo concentrations in this particular model is still 
i 860 under investigation. In contrast to the observations, the CIO 
decreases obtained by B2dM, EMAC, and SOCOL are not 
pole-centered but shifted slightly to lower latitudes. The re- 
maining models (except CAO) produce a very similar CIO' 9 ' 5 
signal. 

was The differences of observed and modeled CIO changes at 
latitudes poleward of 70° N are related to the background 
CIO abundances. Fig 33 shows the zonal mean distribu-1920 
tions of CIO vmrs on 26-27 October prior to the onset of 
the proton forcing. CIO vmrs of more than 0.4 ppbv have 
1900 been observed around 2 hPa in the entire NH with a slight de- 
crease poleward of 70°N. Maximum abundances were found 
at 60-70°N, exactly at the same latitudes where the CIO in^92s 


creases during the following days occurred. Although we 
cannot exclude that the observed CIO in this latitude range is 
affected by gain calibration errors, this coincidence is some- 
how remarkable. In principle, the enhanced CIO abundances 
around 60-70°N can be related to differences in the latitu- 
dinal distributions of daytime OH and O, the first being re- 
sponsible for CIO production and the latter for CIO removal. 

Modeled CIO abundances do not show this enhancement 
around 60-70°N. Simulated CIO vmrs are also generally 
lower by 50% than those observed by MIPAS (except for 
FinROSE and CAO) and exhibit a pronounced decrease to- 
wards the polar night region. In some models (e.g., B2dM 
and EMAC) CTO has disappeared nearly completely at the 
pole. It is thus not surprising, that modeled CIO depletions 
at 70-90°N are less pronounced than in the observations in 
absolute terms. The much stronger modeled decrease of CIO 
towards the polar night region during pre-SPE conditions 
seems to be related to an overestimation of CIO losses. Since 
the sequestering into the CI 2 O 2 dimer is inefficient around 
2 hPa and simulated HOC1 or CIONO 2 distributions before 
the SPE do not indicate a conversion of CIO to these species, 
it is most likely that CIO is more efficiently converted to HC1 
than indicated by the observations. The faster conversion in 
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Fig. 34. Temporal evolution of area-weighterd averages of relative HOC1 changes with respect to 26 October 2003 in MIPAS observations 
and model simulations at 70-90°N, as well as differences between modeled and observed averages. The significance of observed HOC1 
changes (in units of cr) is shown in the upper second panel (from the left). 


the models might be related to the reaction path CIO+OH 
— > HCI+O 2 which has an uncertainty of its rate constant of 
several 100% (Sander et al., 2006). However, also dynami-1950 
cal reasons (i.e., differences in the magnitude of meridional 
1930 mixing) cannot be excluded. 

The temporal evolution of observed and modeled HOC1 
changes at 70-90°N until mid November is shown in Fig. 34. 
HOC1 started to increase rapidly on 29 October, reaching 19 ’ 
values around 0.25 ppbv, and diminished after 1 November 
1935 within a few days. A smaller second increase occurred on 
3 November related to the second, weaker SPE, The simula- 
tions show generally smaller enhancements (approximately 
30% less on average), except FinROSE. This model overes - 1960 
timates significantly the observed enhancements by nearly a 
1940 factor of 2. There, HOC1 abundances remain enhanced after 
the SPE for nearly one week and show a second, even more 
pronounced enhancement around 1 1 November. 

1965 

Figure 35 shows observed and modeled changes of the 
HOC1 zonal mean distribution averaged 29-31 October with 
1945 respect to 26-27 October. The pronounced anti-correlation 
of HOC1 increases and CIO decreases (compare Figures 35 
and 32) indicates that ambient CIO is quickly converted to 


HOC1 via reaction R14 during nighttime in the presence of 
proton forcing. However, HOC1 increases are higher than the 
corresponding CIO losses, resulting in a net increase of ac- 
tive chlorine by approximately 2 ppbv in the observations and 
most of the models. This can be explained by SPE-related 
chlorine activation via reaction R12. FinROSE, however, 
overestimates the chlorine activation by a factor of 3. 

The sharp decline of the HOC1 enhancements after the 
proton forcing observed by MIPAS, and also reproduced by 
most models, must occur in the sunlit atmosphere close to 
the polar night terminator, since losses via reaction R15 are 
negligible in the polar night stratosphere after the SPE. This 
is also the reason for the relatively long lifetime of the HOC1 
enhancements in B2dM where meridional redistribution is 
weak. This is not the case in the FinROSE model. There, 
the long lifetime of the HOC1 enhancements related to the 
SPE, as well as the second buildup around 1 1 November, 
seem to be caused by an underestimation of chemical losses 
ofHOCl. 

CIONO 2 increases, occurring approximately two days af- 
ter the SPE, are attributed to the termolecular reaction 

C10 + N0 2 +M-> CIONO 2 +M. 


(R16) 
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Fig. 35. Zonal mean HOC1 changes during the main proton forcing (29-3 1 October) with respect to 26-27 October in MIPAS observations 
and model simulations. Solid contour lines reflect 0. 1 ppbv steps. 


1970 CIONO 2 is removed mainly by photolysis in the sunlit atmo- 
sphere and, to a lesser extent, by reaction with atomic oxy- 
gen. Due to its pressure dependence, CIONO 2 formation byisss 
reaction R16 is more effective at lower altitudes. Enhanced 
NO 2 availability related to the SPE, however, is increas- 
es ing with altitude, leading to a peak height of the observed 
CIONO 2 enhancements around 3hPa (^36 km). This is 
slightly higher than reported by Lopez-Puertas et al. (2005a)ooo 
who based their analysis on an older MIPAS CIONO 2 data 
version than used here. 

i9so The temporal evolution of observed and modeled 
CIONO 2 changes at 70-90°N until the end of November is 
shown in Fig. 36. The observed enhancements of 0.4 ppbv 
after the SPE remained in the stratosphere for about two 
weeks. After a sudden decrease on 13 November, CIONO 2 
1985 abundances raised again on 19 November, reaching a sec- 
ond, weaker maximum around 22 November. The modeled 
CIONO 2 increases are generally smaller (except CAO, sees™ 
discussion above) and show a different temporal evolution. 
The CIONO 2 underestimation in the simulations, particu- 
1990 larly during the first enhancement starting on 1 November, is 
related to the reduced CIO availability compared to the ob- 
servations. Figure 37 shows the corresponding zonal mean*™ 


distribution of the observed and modeled CIONO 2 enhance- 
ments averaged over 1-5 November. From the observations, 
it is clear that CIONO 2 is principally formed in the polar 
night region where high NO 2 abundances are available and 
no photochemical losses occur. Most of the model simula- 
tions, except CAO, SOCOLi and WACCM, show negligi- 
ble enhancements there. Instead, CIONO 2 formation occurs 
around 70°N, were daytime losses are still small but CIO is 
available, however, with a considerably smaller magnitude 
than observed. SOCOLi and WACCM simulations, which 
have a similar latitudinal distribution of CIONO 2 changes as 
observed, exhibit higher CIO abundances in the polar night 
region than other models. 

The observed temporal evolution of the CIONO 2 changes 
in the second half of November is better captured by mod- 
els based on ECMWF- and MERRA-driven meteorology up 
to the stratosphere (i.e., B3dCTM, FinROSE, KASIMA, and 
WACCM), which hints at a strong impact of vortex dynam- 
ics on the CIONO 2 abundances. Wave-driven vortex ex- 
cursions to illuminated latitudes, alternated by reformation 
of a pole-centered vortex, are mainly responsible for the 
CIONO 2 variability and particularly for the decrease around 
1 6 November. The descending N0 2 layer, formed during the 
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Fig. 36. Temporal evolution of area-weighted averages of relative CIONO 2 changes with respect to 26 October 2003 in MIPAS observations 
and model simulations at 70-90°N, as well as differences between modeled and observed averages. Contour lines reflect 0.5 ppbv steps. The 
significance of observed CIONO 2 changes (in units of a) is shown in the upper second panel (from the left). 


SPE, acts as a reservoir for continuous CIONO 2 formation 
in the following weeks after the SPE. Due to the reduced 
CIO availability in the polar stratosphere towards the end 
of November, additional CIONO 2 buildup is observed only 
around 60-70°N, in agreement with most of the model re-2040 
suits (not shown). B2dM and HAMMON1A, however, show 
very small CIONO 2 increases in the second half of Novem- 
ber. In the first model, this is related to the confinement of 
the NO 2 layer to high latitudes, where no CIO is available. 

In the latter model, strong meridional mixing led to a dilution 
of the SPE-generated NO 2 layer, such that insufficient NO 2 
was available for additional CIONO 2 buildup. 

7 Conclusions 

2050 

We have compared stratospheric and mesospheric composi- 
tion changes observed by MIPAS in the NH during and af- 
ter the Halloween proton event with simulations performed 
with state-of-the-art atmospheric GCM and CTM models. 
The large number of models participating in the intercompar- 20 E 5 
ison exercise allowed for an evaluation of the overall ability 
of atmospheric models to reproduce observed atmospheric 


perturbations generated by SPEs, particularly with respect to 
NO y and ozone changes. This model validation represents 
a mandatory first step towards an accurate implementation 
of particle precipitation effects in long-term climate simu- 
lations. It has also allowed to test and identify deficiencies 
in the chemical schemes, particularly with respect to nitro- 
gen and chlorine chemistry, being relevant for stratospheric 
ozone. 

Observed SPE-related short-time increases of the minor 
species HNO 4 and H 2 O 2 have been identified for the first 
time and are qualitatively reproduced by the simulations. The 
observed enhancements of 0.2 and 0. 1 ppbv, respectively, are 
overestimated by the models on average. Both observations 
and simulations give further evidence for an SPE-induced 
CO depletion. A clear isolation of these chemical losses from 
dynamical variability, however, is difficult. 

In general, atmospheric models are able to reproduce most 
of the observed composition changes. In particular, simu- 
lated SPE-induced ozone losses agree within 5% with the ob- 
servations on average. This excellent agreement is found on a 
short-term scale (HO x -driven) in the mesosphere, as well as 
on a mid-term scale (NO x -driven) in the stratosphere. Sim- 
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Fig. 37. Zonal mean CIONO 2 changes after the main proton forcing (1-5 November) with respect to 26-27 October in MIPAS observations 
and model simulations. Solid contour lines reflect 0. 1 ppbv steps up to 0.6 ppbv and 0.2ppbv steps above. 


ulated NO y enhancements around 1 hPa are on average 30% 
higher than indicated by the observations, while an underes- 
2060 timation of modeled NO y of the same order was found in 
the mesosphere. The systematic behavior in the stratosphere 
suggests that these differences are related to the simulatedoas 
ionization rate profile shape. WACCM simulations without 
inclusion of electron-induced ionization yield much better 
2065 agreement with the observations around l hPa than the than 
the nominal model run including both electrons and protons 
throughout the period of interest. An upper stratospheric ex -2090 
cess NOy production by electron-induced ionization of 5- 
1 0 ppbv could be identified from these simulations, partic- 
2 o?o ularly after the minor second event around 4—5 November. 
The excess ozone loss related to electron- induced ionization 
has been estimated to be around 5%. Again, better agreement 
with the observations was achieved when excluding the elec- 2095 
tron contribution to ionization in the WACCM simulations. 
2075 Our comparisons hence suggest that the modeled electron- 
induced ionization during the Halloween event is consider- 
ably overestimated. This might be related to the uncertainties 
in constraining the electron energy spectrum from available™ 
particle observations. 

The impact of chemical NO losses due to reaction with 


atomic nitrogen (R2) on the SPE-inducedNO y increases has 
been studied in detail. An important dependence of the net 
NOy generation on temperature and background NO x due 
to this mechanism has been identified. In the stratosphere, 
SPE-related NO y increases are reduced (enhanced) by ap- 
proximately 10 % if temperatures were 10 K. lower/higher. 
This behavior might be of relevance for future implications 
of SPE effects on climate when considering a stratospheric 
cooling trend related to climate change. The reduced NO y 
production efficiency related to reaction R2 also implies lim- 
itations for models using family approaches in their chemi- 
cal schemes, since this mechanism of NO y destruction is not 
taken implicitly into account in these models. 

The analysis of the observed and modeled NO y parti- 
tioning in the aftermath of the Halloween SPE has clearly 
demonstrated the need to implement additional ion chemistry 
into the chemical schemes. Short-term HNO 3 increases can 
only be reproduced by model calculations including ion-ion 
recombination between NO 3 ” and H + cluster ions (Verro- 
nen et ah, 2008). The partitioning of HNO 3 and N 2 O 5 in 
the following weeks after the SPE is significantly underes- 
timated by the models that do not include HNO 3 formation 
via water cluster ions (Bohringeret ah, 1983). However, fur- 
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ther work is required to tune the parameterizations of these 
210s mechanisms such that quantitative agreement with the obser- 
vations can be achieved. 

The overestimation of observed H2O2 and HNO4 en-' 80 
hancements by the models hints at an underestimation of the 
OH/HO2 ratio in the upper polar stratosphere during the pro- 
2110 ton forcing. F urther work is required to analyze in detail pos- 
sible reasons for this behavior. The analysis of SPE-induced 2lfi5 
changes of the chlorine species CIO, HOC 1 and C 10 N 0 2 
has shown that the encountered differences between mod- 
els and observations, particularly the underestimation of ob- 
2115 served CIONO2 enhancements, are related to a smaller avail- 
ability of CIO in the polar region already before the SPE. 2170 
In general, the intercomparison has demonstrated that dif- 
ferences in the meteorology and/or initial state of the atmo- 
sphere in the simulations causes an important variability of 
2120 the model results, even on a short timescale of only a few 
days. The model responses to the proton perturbation thus 
show a significant spread. On the other hand, this sensitiv- 
ity of the simulated atmospheric responses to the background 
conditions, indicated by the spread in the model results, also 
2125 implies that the real atmospheres response to proton events™ 
depends strongly on the actual conditions. 

Future HEPPA model-data intercomparison activities will 
focus on the assessment of indirect effects of energetic par- 
ticle precipitation related to polar winter descent of upper 
2130 atmospheric NO x generated by electron precipitation. Thi^ ,as 
is motivated, on one hand, by the higher potential of indi- 
rect effects to influence middle atmospheric composition on 
longer time scale's compared to direct effects (i.e., SPEs) and, 
on the other hand, by its large variability related to dynam- 2190 
2135 ical modulations, making its representation in current atmo- 
spheric models challenging. 
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